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AllSTRACT 

Oxygen steelniaking converter is the single most important unit process in a steel plant and is 
some times called the throbbing hoarl of a steel plant. Introduction of sublance technology in 
the early eighties to analyze carbon and temperature in an oxygen steelmaking converter 
during the blow but approximitlely 3-4 minutes before the end blow proved to be an 
important step in the control of oxygen steelmaking processes. Since the introduction of 
sublance, the accuracy of end poinl prediction (hit rate) at most of the steel plants has 
gradually increased from approxiittalely 60% to 90%, The efforts are now being made to 
narrow down the range of end point carbon and temperature, known as carbon window and 
temperature window. Several models have been reported in literature and these are critically 
reviewed. 

It is generally agreed that tlie rate limiting step for oxidation of carbon 
is mass transfer in metal phase. However, it lias been hard to isolate as to which is rate 
controlling : the mass transfer in bulk metal or the mass transfer in metal droplets ejected 
from jet impact zone. In the present work both tlie possibilities are investigated by developing 
a new model basetl on mass transfer control in which the operational factors like lance 
height, time of blow, coolant additions, slag volume, slag basicity, etc are also incorporated. 
Exhaustive operational data collecleci iVom four different plants are analyzed and a 
comparative evaluation of dilTcrciil models has been done to find out the suitability of 
different types ofmodel.s. d'he average end poinl carbon lies in the range of 0.035-0.06 % and 
average end point temperature lies in the range of 1600-1660 ”C. U is now possible to predict 
end point carbon within a = :t X mill! percent and temperature within o = ± 8 *^0 with the 
help of models developed in the iircsent work. 
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CHAPTER 1 


OXYGEN STEELMAKING PROCESS AND THE 

SUB LANCE 


1.0 Introduction 

Oxygen steel making process accounts for approximately 65% of the steel produced 
in the world. A typical diagram of the physical state of an oxygen steelmaking converter is 
shown in Fig 1. The pure oxygen jet emerging from the lance, placed at center of vessel, 
impinges on metal surface. A stream of metal droplets is ejected due to deflection of the jet in 
the impingement zone. The droplets eventually fall back to metal after their passage through 
slag. The amount and size range of the droplets depend on lance nozzle design and lance 
height whereas the residence time of droplets in slag depends upon the properties of slag. The 
chemical reactions occur primarily in the jet impact zone and in the droplets falling through 
slag. 

In the early stages of its developments, oxygen steelmaking 
process itself was given several names like LD (Linz-Donavitz), EOF (Basic Oxygen 
Furnace), BOP (Basic Oxygen Process), BOS (Basic Oxygen Steelmaking), etc. These names 
were modified when bottom blowing process was invented in which the oxygen gas was 
blown through tuyeres located at the bottom of the converter. The tuyeres needed cooling 
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because of the high temperature generated at the tuyere tip. Depending upon the nature of 
coolant used (viz. fuel oil or natural gas), the new breed of processes were named as Q-BOP, 
LWS, etc. As a further modification of bottom blowing practice, combined top and bottom 
blowing processes were invented and today more than 30 commercial names exist for the 
variety of oxygen steelmaking processes. For the sake of simplicity, however, they all are 
now grouped under a single name of ‘Combined blown processes’. 

The specific combination of top and bottom blowing (i.e. the flowrate 
and nature of gas through top and/or bottom) adopted at a particular plant depends upon the 
starting composition of hot metal, capacity of converter, quality of lime, reffactoiy life, 
desired end point composition of steel, economies of operation, etc. As a result, the 
steelmaking process at each steel plant has its own unique features and the control models 
developed at one plant are rarely applicable to the situation at another steel plant. In a more 
recent development (mid eighties), a separate lance called “Sublance“ is introduced into the 
converter to measure temperature and take a sample of metal, approximately 3-4 minutes 
before the end of the blow. Corrective actions are then taken based on analysis of sample 
obtained by sublance measurement and this has significantly improved the accuracy of 
control of obtaining desired end point temperature and carbon of liquid steel. The fine tuning 
models so developed are called “sublance models”. 

1.1 Scope of Present Work 

The objective of the present work is to develop sublance models (as described 
above) for the prediction of end point temperature and carbon of liquid steel. Different 
models available in literature are reviewed in Chapter 2 and tested on operational data in 
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Chapter 3. A new model is also developed in the present work based on the assumption that 
decarburization primarily occurs in the recirculation zone of jet impingement. Data have been 
obtained from different steel plants (Chapter 3) to compare if the same nature of variables can 
be used to characterize the last stages of oxygen steelmaking process. The contribution of 
droplets to decarburization vis a vis the dec^urization in recirculation zone and the 
performance of different models is also discussed in Chapter 3. The main conclusions of the 
investigation are summarized in Chapter 4 along with suggestions for further work. 
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CHAPTER 2 

REVIEW OF MODELS FOR END POINT CARBON AND 
TEMPERATURE PREDICTION USING SUBLANCE DATA 

2.0 Introduction 

One of the key objectives during end blow period (i.e. the period after sublance 
measurement) is the control of decarburization rate. A lot of research has been done in this 
area and several models, both static and dynamic, have been proposed which are suited to 
specific steel plant conditions. A brief summary of the fundamental basis of various models 
based on sublance measurements is presented in section 2. 1 to section 2.5. 

2.1 Model by Lee and Rao 

Lee and Rao [1] have developed a mixed control model for decarburization in 
levitated droplets assuming gas phase mass transfer and dissociative adsorption of CO 2 (or 
interfacial chemical reaction) by studying the kinetics of decarburization in high carbon 
regime (0.91 1 - 2.479 wt % C) at 1973 K. The rate equation assuming mass transfer control 
in the gas phase as 
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1200.A 

pV 


RT, 


•ln(l + P*co 2 ) 


( 2 . 1 ) 


where, 

Cc = wt % of carbon in the melt 
A = interfacial area (cm^) 
p = density of the melt (g/cm^) 

V = Volume of the melt (cm^) 

kg = average mass transfer coefficient in gas film (cm/sec) 
Tf = average gas film temperature ( °K) 

P'’co 2 = bulk pressure of CO 2 . (atm) 

If interfacial chemical reaction control is assumed as a rate controlling step, then 



dt 


1200 ^ 

pV 


l^KejCs 

1 +KfC ^ 


( 2 . 2 ) 


where, 


sulfur 


metal. 


kR = rate constant (mol/mVsec/atm) 

K = adsorption coefficient of Sulphur 

0u = the jfraction of the surface sites not available for the chemisorption of 

f = activity coefficient of sulphur in the surface layer and in bulk 

Cs = bulk concentration of sulphur (wt %) 



P‘co 2 = interfacial pressure of CO 2 . 

On assuming mixed control of mass transfer and interfacial chemical reaction, the 
decarburization rate is 



1200 ^ k^k^iy + KOJCs) 

pV '(k\+k^) + {k\+k R9^)KfC, 


ln(l + P^o,) 


(2.3) 


where, 

* RI^ 

The experimental results of Lee and Rao showed that 

• Liquid phase mass transfer is not a rate controlling step in the high carbon regime. The 
initial carbon concentration has no discernible effect on the decarburization kinetics. 

• Sulphur has a clear and reproducible retarding effect on decarburization. This effect is 
most pronounced in the range of 0 to 0.05 mass % sulphur. 

• The Pco 2 of the gas mixture has a marked effect on the decarburization kinetics. 

• The flow rate of the gas mixture has a small but finite effect on the rate of 
decarburization. 

In low carbon range (C< 0.4 mass %) equation assuming liquid phase mass 


transfer as control step is 



where, 
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Rc = rate of carbon transfer (mol / sec) 

kc = mass transfer coefficient of carbon (m/sec) 

nc = concentratin of carbon (b=bulk, i=interface ) 

Distin, Hallett and Richardson [2] have assumed in their study on the decarburization of 
liquid Fe-C alloys with O 2 and CO - CO 2 mixtures that the concentration of carbon at the 
surface becomes virtually zero at the instant of oxide formation. When the oxide appears on 
the surface Eq. (2.4) is reduced to 




k^A j 
— n c . 
V 


Or 


-4%C] 

dt 



(2.5) 


where, 

A = Surface area of droplet (m^) 
V = Volume of droplet (m^) 

For spherical droplets [1] 



( 2 . 6 ) 


where, 

D = difflisivity of carbon 
a = radius of the droplet 

Lee and Rao calculated the kc for their experimental conditions as follows: 


T = 1973K 
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& = 0 . 1 % 


a = 0.2875 cm 
D = 3.24X10'^ cmVsec 


k - 5^3.24x10"^ 

a ~ 0.2875 


56.348;cl0 ^cm/sec 


= 56.348 X 10'^ m/sec 

They found that mass transfer in the liquid phase is enhanced by about ten 
times due to electromagnetic levitation as compared to that at stationary condition. This 
enhancement is brought about as a consequence of motion in the liquid phase. A situation, 
similar to those for levitated droplets, may be assumed to prevail during decarburization in 
converter where droplets are ejected from jet impact zone at high velocity. 


2.2 Model by Chou et al. 

Chou, Pal and Reddy [3] have developed a general model for the 
decarburization process. If it is assumed that hot metal at 1400°C with a high carbon content 
is initially charged (for instance, liquid iron containing 4% carbon) in the converter with 
some scrap, then for the sake of simplicity bath temperature may be approximated to vary 
linearly with carbon content as follows: 

T= 1873 -50[%C] (2.7) 

where, 

T = the absolute temperature of the bath 

[%C] = carbon concentration in the both expressed in wt. pet. 

In an actual converter, bath temperature can be predicted at each time step by 
dynamic heat and mass balance. When oxygen blows into the liquid metal, the carbon in hot 
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metal will be oxidised by oxygen to form carbon dioxide. At an elevated temperature, the 
CO 2 will react again with bath carbon according to the Boudouvard reaction 

C02+[C] =2CO (2.8) 

where, 


[C] = the carbon dissolved in the hot metal. 

In accordance with the mass action law, we have 




CO 

Pco, 


c 


(2.9) 


where, 

Pco , Pcch = the partial pressure of CO and CO 2 respectively, 
ac = activity of carbon in hot metal. 

K, the Equilibrium constant, is a function of temperature and can be derived 
in terms of standard Gibbs free energy data. 


K = exp [15.30 -(16759/r)] 


( 2 . 10 ) 


Substituting Eq (2.7) into Eq (2.10) and combining with Eq (2.9) yield 


exp 


15.30- 


16759 


1873-50[%C]J 




CO 


PcOi -^c 


( 2 . 11 ) 


If the oxygen blowing into the converter bath is not pure but contains some inert gases such 
as nitrogen, argon etc. and the total pressure is assumed to be P, we would have 
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Pco ~ P- Xco 

(2.12) 

Pc02 = P- Xc02 

(2.13) 

Pinr P. X inr 

(2.14) 


Where xco, xco 2 , xinr represent the mole fractions of CO, CO 2 and inert gas in the 
decarburization zone, respectively, and 


Xco + XC02 + Xinr = 1 


(2.15) 


On the other hand, the activity of carbon can be expressed in terms of the 
products of [%C] and activity coefficient, fc i.e. 


ac=[%C]fc (2.16) 

Substituting Eqs (2.12) to (2.16) into Eq (2.11) generate^ 


exp 


15.30- 


16759 

1873-50[%C] 


Px^co 


(2.17) 


Let two other quantities be defined as Xco and Xco 2 (which are relative 
contents of CO and CO 2 in CO + CO 2 ) 


Xco ~ Xco/ (XCO+XCO 2 ) ~ Xco/ (1 - Xinr) (2-18) 

XC02~XC02/(XCO+XC02)='XC02/(1 -Xinr) (2-19) 


Substituting Eqs (2.18) and (2.19) into Eq (2.17) we obtain 
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exp 


15.30- 


16759 

1873 - 50 [%q 


PX^co(\-Xj.^) 

(l-Z^J[%C]./e 


( 2 . 20 ) 


The above relation gives the dependence of relative content of the carbon 
monoxide on the carbon concentration in the hot metal during the BOP operation. Based on 
this solution the relative mole fraction of CO can be calculated if carbon content and fc are 
known. 


Bensel et al [4] have used the following expression for fc. 


fc=a[%C] (2.21) 

where, 

a = 0.167 when [%C]>1% 

a = 0.2 when [% C] < 1% 

These are approximation formulae and they completely neglect the effect of temperature 
on fc. Rist and Chipman [5] have given a Equation, which explains the effect of temperature 
on activity coefficient of carbon 


log fc = 0.1666 [% C] - 0.01585 [% Cf + 9.9613 X lO*'" [% Cf (T-273) 

+ 3.0246X10-^ [%Cf (T- 273) (2.22) 


The decarburization rate depends on the oxygen blow rate into the furnace. If 
the rate of oxygen blow is dno 2 /dt (mol/ min) and the decarburization rate dnc/dt, then from 


mass balance considerations 
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= (2.23) 

Where the factor (1 - Xinr) reflects the purity of oxygen. The negative sign stems from the 
decreasing carbon content during decarburization. It is implicitly assumed that oxidation of 
other elements like P, Mn, Fe, S is negligible in high carbon regime and this is approximately 
so. 

Expressing in terms of [%C] and the flow rate Qo (Nm^/ton) the above 
equation becomes 


4%C] _ 

dt 1.87xl0-'(l-0.5^co)i^ 


(2.24) 


where, 

Wt = bath weight in gms 
[%C] = weight percent of carbon 

When carbon content is less then 0.3 , the decarburization will be controlled 
by mass transfer of carbon in the melt instead of oxygen and it can be expressed as 

dnc/dt = - kd. A. (C*’ - C) (2.25) 

where, 

C’’ , C = bulk and interfacial concentration of carbon, mol/cc 
kd = mass transfer coefficient (cm/sec) 

A = bath / oxygen interface area (cm^) 

Since the reaction between carbon and oxygen at high temperature is very fast, Ci is 
legligible (except towards end of the blow) and hence Eq (2.25) can be written as: 
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dnc/ dt = -kd.A. (C^) 


Or 


4%C] 

dt 



where. 


V = melt volume 


(2.26) 

(2.27) 


2.3 Exponential Model 


2.3.1 Model by Bessho et al. 

Bessho et al.[6] have developed a dynamic model, based on exponential 
dependence of decarburization rate on carbon content . This model has been applied to 
estimate carbon content and temperature of Q -BOP heats after sublance measurement. 
Exponential model is expressed as 


dC 

dO^ 


= ^ + z'exp(^) 


(2.28) 


where . 


T=-^ exp| 


C -C 
'-0 


and (j> 


-1 


C, - Co 


Then Eq (2.28) becomes 
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dC 

dO^ 


= & l-exj 


( c-c 




c -C 


oy 


(2.29) 


where. 


^ = theoretical maximum decarburization rate 

Co = Critical carbon content below which decarburization does not substantially occur 
and is found to be 0.01% for Q- BOP. 

Cp = a parameter derived statistically based on operational variables. 

Value of can be calculated for a converter with oxygen flow rate Qo = 600 
Nm^/min and weight of bath Wt = 300 tons as follows: 


1 = 


600 

22.4' 


-x\2xl 


300x1000 


= 0.214%C/min 


By integrating the Equation (2.29) 


I 


-dC 


1-exp 


Q-C 

Q-Coy 




Or 


J C„-C ~ ^^2 




Cq-C 

Let t = 


log^t = 


Cq-C 

Cp-Co 


By differentiating, we get 


dt d C, C 

t ■ Jc[c^-Co“C^-Co 
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Or dC = 

t 


Substituting for dC value, we get 

f 1 -(C,-C„) , 

-#Ao, =JYr7- — ^ 


Or, $iO,=(C,-Q,)l[y^+y]r4 


= (C -C.).(-ln[l-(] + ln[r)) 


( C. -C^ ^ SizSS 

(C, - cA -Inl 1 - exp-^ 1 + ln| 

|_ V ^oj \ J_ 


[ r -c ] r -c 
(C,-Co)-ln 1-exp-^^ 


Taking the integration in the limits of Csl to Cfa, we get 


^.AO,=(C -Co) -ln 


r -c '\ 

C,-C« . Q-C,, C,-C„, 


1-exp 


C.-Crr C,-C, C,-C. 

c,-c. 


f c ^ 

‘ “’’c^-c, Q-C,.. 

|.AO, = (C, - CJ -In. ,/_c +-F^ 

V S ~ ^0 y 
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AOj = 


(Q-CJ 



( C -C 

1-exp ^ 


In 

C -C 

'-0 

1 C/a ^SL 

C -C 
l-exp^® 

c -C 

Cp-C, 

_ 

\ P ^ J 



(2.30) 


where, 

CsL = carbon content determined by sensor lance 
Cfa = the aimed carbon content 

A O 2 = the amount of oxygen required (per ton of charge) to reduce the carbon content 
from CsL to Cfa . 

The carbon content at any arbitrary time after the measurement by the sensor lance 
can be calculated as follows. 

From the Eq(2.32) , we can get 


Q-C, 


^fa ^SL 


In 


1-exp 


^0 ^fa ^ 

C,-Co 


1 Cq-Csl 
1-exp— 2 ^ 

C -C 

'-'p '-0 


Or, exp[ 


■^.AO, C,,-C 


1-exp 


SL 


C -C C -C 
'^0 ^0 


]= 




1-exp 


-C ^ 

*^0 '^SL 
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where 


R = exp 


V 0 ;,-Co J 
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where AO 2 (Nm^ /ton of charge ) is the amount of oxygen blown after the measurement by 
sensor lance. The value of Cp in the above equations is estimated from multiple linear 
regression (MLR) analysis which involves the carbon content and temperature of the bath 
monitored by the sensor lance.The MLR equation to estimate Cp can be expressed as 



m 

i=l 


(2.32) 


Thus, 

ao,a i = regression constant and regression coefficients. 

Xi = operational variables as shown below 
The operational variables used by Bessho et al. For regression analysis are 
Xi meaning 

1 real carbon content at blowing end 

2 real bath temperature at blowing end 

3 hot ratio (=wt of bath /wt. of scrap ) 

4 carbon content by sensor lance 

5 bath temperature by sensor lance 

6 amount of lime used after one point check by sensor lance. 

7 total amount of lime and dolomitic lime used 

8 calculated basicity 

An attempt has been made in Chapter (3) to apply the above relationship to the present work 
data (301 heats). It may however be noted that it is not correct to include the operational 


variables like 
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Xi = real carbon content at blowing end 
X 2 = real bath temperature at blowing end 
X 4 = carbon content by sensor lance 

as independent variables in the regression Equation to estimate the dependent variable ‘Cp’. 
This is because Xi, X 2 , are not known a priori at the time of sublance measurement and X 4 
is already used in the Eq (2.30) to calculate Cp 

Bessho et al. used a semi empirical equation to estimate the temperature rise. 


where. 








(2.33) 


Tf , Cf = the temperature and carbon content of the bath at blowing end respectively. 
Tsl ,Csl = the temperature and carbon content of the bath measured by the sensor 
lance 


AO 2 = amount of blowing after sensing ( NmVton of charge) 

P, Q, R = Parameters calculated from multiple regression analysis by using actual 
operational data such as lance height, slag volume , etc. 


2.3.2 Model by Byun et al. 

Byun et al. [7] employed a simplified form of exponential model Eq. 

(2.29), as: 

= (2-34) 


Wst = Weight of molten steel (kg). 
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Here the term P is same as l/(Cp -Co) and a is same as ^ in Eq (2.29). It is implied that p 
represents combined effect of Cp and Co on decarburization. 

Integrating Eq (2.34) from sublance carbon content Cs to aimed carbon Ce gives the oxygen 
consumption from the in - blow sublance measurement to the blow end 



c. 










-W 

d[C] = 


a 


/- 


1 


-yS(C-Co) 


.d[C] 


Let t = 

In (t) =-p (C-Co) 

By differentiating, we get 


(2.35) 


dt 




dC = 


-1 dt 

T~ 


On substituting dC value 


-W r I -I dt 

AO. 

a Jl-t P t 


Or, 


AO. 


-W^ 7 1 dt 


2,caZ 


= Z1L^ f- 
cc.p ^ 1 


■t t 


a.p It 1 - 




^.[in(o-Ki-o] 


By substituting t value, we get 
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= - ln(I -e-*'-'-*)] 

a.p c. 

W 

Or NO = ^ 

VJl, i^'^2.cal o' 

a.p 


AQ-C,)-ln 




\-e 


-PiC.-C,) 


JJ 


(2.36) 


The effect of other operational variables on oxygen consumption can be incorporated through 
multiple linear regression as 


AO. 


2. cal 


w 

" St 

a.p 


AQ-Q-ln 






(2.37) 


where, 

Yi = regression coefficien 
Yo = regression Constant 
Xi = operational variables , 

= Standard operational variables 

2.3.2. 1 Temperature calculation model 

the heats evolved in 

(2.38) 

(2.39) 


Temperature at the end of blow can be calculated by considering 
various reactions 

[C] + l/202(g) = CO(g) , AHi 
CO(g) +1/2 02 (g) = CO 2 (g) , AH 2 
Fe + 1 / 2 O 2 (g) = (FeO) j AH 3 


(2.40) 
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ViOi = [0] , AH 4 (2.41) 

Where AHi is the heat produced from each reaction (Kcal / Nm - O 2 ) 

Over all heat balance (excluding heat losses) can be written as 

Cp Wst dT = AH, dOc + AH 2 dOco + AH 3 dOpe + AH 4 dOd^ (2.42) 

where, 

Cp = Specific heat capacity of molten steel (Kcal / C) 

Oc = oxygen combined with carbon to form CO (Nm’ - O 2 ) 

Oco = oxygen combined with CO to form CO 2 (Nm^ - O 2 ) 

Ope = oxygen combined with Fe to form FeO (Nm^- 02 ) 

Odis = oxygen dissolved in liquid steel (Nm'^- 02 ) 

In Eq(2.40), the heat of Fe oxidation and the heat produced by oxidation of CO to CO 2 
are not greatly different at steel making temperature, i.e. AH 2 » AH 3 . In the later stage of 
EOF blowing process , it is generally accepted that the product of dissolved carbon and 
oxygen in liquid steel is constant. 

[C] x [0] = K '4 (2.43) 

Therefore the amount of dissolved oxygen can be expressed as a function of 
carbon content in steel . 


dO^x 


22.4 
32 100 


22.4 W, K\ 
32 100 [C]- 


4C] = 


(2.44) 


Assuming that the carbon removal from steel melt mainly occurs by Eq(2.38), the oxygen 
consumption for decarburization can be obtained as 


dO^ = 


12 100 


d[Q = -K,d[a 


(2.45) 
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If an assumption is made that all lance oxygen is consumed by reactions (2.38) to (2.41) , 
d 02 = dOc + dOco + dOdis, 
the term dOpe + dOco is expressed as follows. 

dOpe + d Oco =d 02 - (dOc + dOau ) = d 02 + Ki d[C] + K 4 d[C]/[C]^ (2.46) 
If Hi are assumed to be independent of temperature , the following differential equation 
can be obtained from Eq(2.42) through (2.46). 

dT = ( 5i + 62 / [Cf ) d[C] + 53d02 (2.47) 

where, 

51 = ( AH 2 - A Hi )Ki / ( Cp . Wst ) 

52 = (AH2-AH4)K4 / (Cp.W^) 

63 = AH 2 / (Cp.W«) 

Integration of Eq(2.47) from in blow measurement to the end blow gives a linear 
function of temperature rise. 

T. c, C, Oj,. 

AT^ = ^dT = S,jdC + S,j:^+S, jdo, 

= 5,(C,-C,) + 5,(^-^) + 5,AO, (2.48) 

Where Ts , Te are temperature measured at in blow and blow-end and O 2 is 
total oxygen blow between these measurements. 

In order to fit this model better to the actual temperature data, parameters 6 i in 
Eq(2.48) are modified by the multiple regression analysis as discussed in the decarburization 
model. In the present work operational variables other than the amount of oxygen and carbon 
content were also included in this model such as the steel temperature at sublance 



measurement, slag volume, HMR, lining life, lance height and the amount of iron ore used. 
Therefore, the equation for the temperature increase can be rewritten as: 


at: 


actual 


= (C, - C J + 5 , (^— ^) + ^3 AO, + 2 z, ^ + ;}ro 


(2.49) 


where, 

Xi , Si = regression coefficients 
Xo = constant 
Zi = Operational variables 
Z,. = Standard operational variables (mean values ) 

2.3.2.2 Feedback model 

There will always be deviation of oxygen quantity and temperature rise for a heat 
•which are not accounted for by the explanatory relationships such as Eq(2.37),(2.48) and -thus 
some part of those deviations will remain unpredictable. In order to minimise the error 
various steps can be taken. The model can be modified to include the perturbing factor which 
may be explained by the selected model parameter. The model parameter for upcoming heat 
can be adjusted by a time series analysis of the posterior calculation of actual turndown data 
of preceding heats. In other words the new model parameters for upcoming heat can be 
expressed as a function of old parameter of preceding heats. 

One of the most important concerns in developing feedback model is how to 
select the 'adapted model parameters. In this work, various model parameters, 6i for 
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decarburization model and (j)j for temperature calculation model were employed and the 
general form of equations used are 


Aa = e, 


w 1 -t-"' ... 1 

1 


(2.50) 


Ar = ^i, 


(C - ^ Aa + X 

^ t' ^ V / 


(2.51) 


In addition to above, the following time - series forecasting methods can be tested to obtain 
appropriate feed back model. 


Method 1 : Single moving average 

This technique consists of taking a set of observed values, finding the average of 
those values, then using that average as the forecast for the next period. 

Ft., =(X,+ X,.i+ + X..N.1) /N (2.52) 

where, 

F t + 1 = the forecast value for the next period 
Xt = the most recent observation 
N = the number of historical data. 

Method 2 : Exponential smoothing method. 

This forecasting is based on weighting the most recent observation with a 
weight of value, K* whose value is between 0 and 1 and weighting the most recent forecast 
vith a weight of value, (1 - K*). 
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Ft+i =K*X, + (l.K*)Ft (2.53) 

where, 

K* = Smoothing factor 

Method 3 : ARMA method 

In this model, the series to be forecast is expressed as a function of both previous values 
of the series ( anti regressive term ) and previous error values from forecasting ( moving 
average term ) . ARMA (2, 1) model used in Byun et al. model was 
Ft+i = A-i Xt + ^2 Xt-i - p-i et (2-54) 

where, 

A,i ,|j.i = coefficients. 

et = the most recent error from forecasting . 

Of the above three methods, the exponential smoothing method has been reported to give 
the best result with a K* value of 0. 1 . 

The procedure for On-line application of the model is as follows: 

The in- blow sublance measurement is made when 80% of total oxygen amount, 
presumed by the static model has been blown. At this point, the decarburization model starts 
a calculation of the oxygen amount to be blown for the aimed carbon content. The 
temperature increment in this dynamic period is then calculated based on the amount of 
oxygen to be blown. If the calculated end- point temperature is higher than the aimed one, the 
amount of coolant to be added is determined. If the calculated temperature is lower than the 
aimed one, the heat has to be over blown. The bath temperature and the carbon content in 
steel at turndown are used for adapting the model parameters for the next heat. 
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2.4 Model by Szegedi et al. 

Szegedi et al [8] have developed a model for the adjustment of carbon content of 
steel by dividing the decarburization process into 3 periods, namely 

I . Thermodynamic period 
n . Reaction kinetic period 
III . Chemical period 



The formulae to calculate the decarburization rate and time taken for the blow for each period 
are as follows. 

The decarburization rate for the period I is obtained on the basis of operational data as 

Vc,i = 0.1071.1 [-0.1421 .Cb - 0.7366.SiB - 0.0425.i + 1.4987 ] (2.55) 


where, 


Vc,i = rate of carbon oxidation(%C/min) 



1 


= intensity of blowing (m^ 02 /t steel, min) 
= carbon content of charge (%) 

= silicon content of charge (%) 
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Cb 
Sie 

Time taken for computing the period I is given by 


1 1 = 100 . SiB [i.(2.03 1 .Cb + T.Sis + 0. 106 .i - 6)]min 


Similarly for the reaction kinetic period n 

Vc,n= 0.1071. i 

and time necessary for second period is 

tn =9.375 (Cb-0.046 ) / (-1.942 - Si.A/(i . B)) 

where, 

A = Cb + 5. 1728. SiB + 0.2981. i -10.52426 
B = 0.4217 - 0. 1426. Cb - 0.4918 - SIb -0.00744.i 
The third period is chemical period for which the rate of carbon oxidation is 


Vc,in = 12.89[%C] . exp (-6590.73^') 


The regression analysis of operational data has shown that 


7323 

2 1422 - r%Cl 0208./ + 0.046 . 

= 128 ^ [%C] 


(2.56) 


(2.57) 


(2.58) 


(2.59) 


1289 


(2.60) 
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The activation energy of the reaction is 

= -R.6590.73 J/mole. (2.61) 

Total time for blowing t = ti + tu + 1 m 

By equating the decarburization rates periods II and III (i.e. VCu = VC m ) they have derived 
a formula to calculate the critical carbon constant [Ck] . 

Ck = 0.208.(i) + 0.046 . (2.62) 

They have experimented this model in the range of 0.01 -4.0 C wt % .and have 
shown that the adjustment and the more exact control of carbon content ensure a significant 
saving of time and material costs. 

2.5 Models Based on Exhaust Gas Composition 

2.5.1 Model by Fukumi et al 

Fukumi et al [9] have developed a refining control model based on exhaust 
gas information measured by mass spectrometer. This model consists of calculation models 
that use sublance information to estimate various values as follows. 

• Calculation and estimation of the oxidation degree ( WO 2 ) 

• Estimation of (T.Fe) 

• Estimation of decarburization 
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• Estimation of temperate in steel 

• Estimation of oxygen content in steel. 

2.5. 1.1 Calculation and estimation of oxidation degree (WO 2 ) 

WO2 is defined as Eq( 2 . 63 ) based on oxygen balance calculated by refining conditions 
(top and bottom gas flow rate, amount of flux, etc.) and exhaust gas information (analysis, 
flow rate) 

WO2 = [ J (fi(Q,F) - fo (V,A))dt ] / WM ( 2 . 63 ) 

where, 

fi = input oxygen 
fo = output oxygen 
WO2 = oxidation degree 
Q = top blow oxygen 
F = amount of coolant 
V = 02,C0, CO2 volume in exhaust gas 
A = sucked air or slopping gas volume. 

,W02 indicates oxygen in slag and steel, and this value is the most important parameter in this 
model. The volume of air sucked fi'om the converter mouth and slopping gas volume are 
evaluated with reference to bottom blowing gas, balance of air, Ar and N2. Analysed values 
of exhaust gas have a time lag behind the actual reactions because of gas transport time. At 
the final stage of blowing when WO2 is rapidly increased, this time lag must not be neglected. 
Thus, reacted volume can be calculated at real time by the estimation of WO2 fi'om Eq( 2 . 64 ) 


which is calculated fi'om actual data. 
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dW02/d02 = a { 1-exp [WO2/3]} 


(2.64) 


where, 

WO 2 = Oxidation 
O 2 = Input oxygen 
a = Constant 
p = Parameter 

2.5. 1.2 Estimation of total iron 

Since WO2 indicates oxygen in slag and steel , WO2 has a close relationship 
with total iron (T.Fe) in slag. In other words, the change of (T.Fe) can be known from the 
change of WO2. Therefore, from (T.Fe) at the first sublance and the increasing amount of 
WO2 after the first sublance, the value of (T. Fe) at the final stage can be estimated in real 
time with reference to Eq. (2.65) - (2.67) 

(T . Fe)s = f (QB, FO 2 , Cs) (2.65) 

Ws(t) = a W02(t) + Wso / { 1 - p (T . Fe)s / 100 } (2.66) 

(T.Fe) (t) = (T. Fe)s + y WO 2 (t) / Ws ~t (2.67) 

where, 

(T.Fe)s = (T. Fe) at in - blow sublance 
QB = bottom blowing gas flow rate 
Cs = [C] in steel at in - blow sublance 
WO 2 = Oxidation degree 

Wso = Slag volume except (FeO), (MnO) , (P2O5) 
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Ws = Slag volume 
t = time 

a , P , y = Parameter 


2.5.1.3 Estimation and prediction of the amount of decarburization 

The amount of decarburization utilised for the estimation of steel temperature 
and prediction of slopping is calculated by Eq (2.68) based on carbon balance from refining 
conditions and exhaust gas information. 

AC = [ J / (ho(V . A) - hi(F)dt] / WM (2.68) 


where, 

ho = Output carbon 
hi = Input carbon 
AC = Amount of decarburization 
V = CO, CO 2 volume in exhaust gas 
A = Sucked air volume or slopping gas volume 
F = Amount of coolant and flux 

At the final stage of blowing , the estimation of decarburization , amount that 


contains time lag is calculated according to the efficiency of oxygen for decarburization. 
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2.5.1.4 Estimation of steel temperature 

The heating of the steel at the final stage of blowing has a close relationship 
with the generated amount of (T.Fe); namely the increased amount of WO 2 and amount of 
decarburization. 

Therefore, heating after sublance is estimated by Eq. (2.69). 

dT/dOz = a . dW02 / d02 - p . dC / d02 - y WsWM + 5 (2.69) 

where, 

T = Temperature in steel 
O 2 = Input oxygen (gas + coolant) 

WO 2 = Oxidation degree 
C = Carbon content in steel 
Ws = Slag volume except (FeO) , (MnO) , (P2O5) 

WM = Weight of steel 
a , 3 , y , 6 = Parameter 

2.51.5 Estimation of dissolved oxygen content in steel 

Oxygen content in steel is estimated fi-om the oxygen balance Eq (2.70) 
calculated fi-om (T.Fe) which is evaluated fi-om exhaust gas information and the steel 


temperature. 

log {[0]} = log { a (T. Fe)} -h b/T Y 


(2.70) 
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where, 

[O] = Oxygen in steel 
T = Temperature in steel 
a ,P ,7 = Parameter. 

2.5.2 Model by Iwamura et al 

Iwa mura et al.[ 10 ] have developed a new endpoint control system of EOF with an 
auto-tuning function for the model parameter based on the measured exhaust gas 
composition. This system can accurately estimate the steel temperature and the carbon 
content at the blow-end through one campaign without changing of the model parameter. The 
dynamic control model consist of an oxygen balance equation and heat balance equation. The 
amount of oxygen required in a heat is expressed by the sum of the amount of oxygen based 
on the change of carbon. 

Fo, _c^)+a, In^+a, +Zl,(_X, - X,)+FB^ (2.71) 

"CT '-fiP 


where, 

ao,ai = coefficient of oxygen balance 
CsL, Cep = carbon contents at sublance measurement and endpoint 
Fo 2 = the amount of oxygen after sublance measurement 
11= coefficient of oxygen in fluxes 


Wsub = amount of flux 
O2 = the amount of oxygen flxues 



35 


WsT = weight of molten steel 
Xi = factor i of operating condition 
X^ = mean value of Xi 

li = coefficient of statistical model 
FBc = adjustment of oxygen balance from the feedback 

In the heat balance model, the fundamental equation has been changed from 
the carbon content to the oxygen basis . In situations where the carbon content of the metal at 
the endpoint has little variation, temperature variation has stronger relation with the amount 
of oxygen than the change of carbon content. The heat balance is represented by the 
following equation. 

Tip -2k = -1,)+FA- (2.72) 


where, 

Tsl,Tep = temperature at sublance measurement and at the endpoint . 
bo, b2 = coefficient of heat balance . 

Y = coefficient of cooling capacity 
li = coefficient of statistical model 
Tsub = cooling capacity of fluxes 
Fbt ~ adjustment of heat balance from the feedback 


2.5.2. 1 Auto-Parameter -Tuning function 
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Usually the decarburization efficiency is changed during blowing for each 
heat. This fact shows that in the oxygen balance equation with a fixed parameter it is difficult 
to obtain a highly accurate estimate of steel carbon content and temperature at endpoint 
through one campaign . To solve this problem, Iwamura et al. have added an auto-parameter 
-tuning function of parameters in oxygen balance Eq.(2.71) using exhaust gas information 
from a mass spectrometer. The parameters of ao and a; are adjusted according to the following 
Equations. 

First of all, Eq(2.73) is altered to the following Equation . 




dC 


1 




(2.73) 


The decarburization efficiency is expressed with the parameter of ao and a; and the estimated 
steel carbon content . On the other hand, the actual decarburization efficiency is calculated by 


(X = Cj.” 


Ao. 


(2.74) 


where , 

G = gain 

a = decarburization efficiency 

i = sampling number fi'om sublance measuring time to blow-end . 

ACi is the change of steel carbon content and AOi is the amount of oxygen during the 
sampling period At as shown below. . 
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AQ =Ci-Ci.i 


(2.75) 


Ao, = 


t’+^ 

jFo,(r)Jr 


r' 


(2.76) 


The amount of decarburization in the exhaust gas during At is calculated from the following 
expression using the exhaust gas information. 


AWr 


= y [C0{T) + C0.{'^)}.0{r).\2 ^ 

~ •* "j/cnnoo A 1 nnn '' 


3600.22.4.1000 


(2.77) 


From Eq. (2.74) and .(2.77) we can obtain 


AfT... 


a = G.— ^.100 

Aa 


(2.78) 


when C is infinity, the term of 1/C in the Eq.(2.73) can be neglected and the parameter of ao 
is expressed from Eq.(2.73) and (2.74) can be shown as 



(2.79) 


In on-line operation, the parameter of ao is decided by the average of the actual 
decarburization efficiency during the maximum decarburization period before sublance 
measuring time as shown in the following equation 
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1 


(2.80) 


where , 

ct max = decarburization efficiency in maximum decarburization period 
m = constant of sampling number 

On real time after sublance measurement, the parameter of ai in Eq.(2.73) is calculated as. 


^ = 


a, 


(2.81) 


The parameter of aj changes in proportion to the decrease of the decarburization 
efficiency and this function can accurately estimate the steel carbon content. 

2.5.2.2 Feedback calculation model 


The feed back calculation model has been developed in order to adjust the 
errors in the dynamic control model. The amount of feedback is calculated by considering the 
trend in the errors of the dynamic control model in past some heats measured by sublance at 
blowing and endpoint. The rule of feedback is classified with three cases by the trend of the 


error as shown below. 
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■ . 




Case 

The trend of the error of the model 

Th amount of feedback 

1 

Mdlerro(i- 1 ) . Mdl erro(i)>0 

FB(i) = FB(i- l)+GFB(erro(i)+erro(i- 1 ))/2 

2 

Mdlerro(i-l).Mdlerro(i)<0 and 

|erro(i)|-|erro(i-l)>H 

FB(i) = FB(i-l)+GFB.erro(i) 

3 

Mdlerro(i-l).Mdlerro(i<0) and 

|erro(i)|-|erro(i- 1)1^ 

FB(i) = FB(i-l) 


where, 

FB(i) = feedback value 
Mdlerro(i) = error of the model 

erro(i) = eror of the feedback model 
Gfb = feedback gain 
H = on-detective value 


In the oxygen balance equation, the feedback value is adjusted according to 
the condition of the sublance measured carbon content at end point in addition to the above 
rule of the feedback. This is because the sensitivity of the carbon content for the amount of 
oxygen varies with the sublance measured carbon content in the blowing. 

The relationship between carbon content and amount of oxygen can be defined in the 
following equation 


^EP ~ -^^8 


\WsrJ 


+J, 


where. 


(2.82) 


Si,S2 = coefficient 
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The sensitivity of the carbon content to the amount of oxygen is related to the following 
equation. 



( 2 . 83 ) 


where, 

8 = the sensitivity 

The actual feedback value is adjusted with the sensitivity as shown in the following 
equation. 


FB’(i) =h . FB(i) 


( 2 . 84 ) 


where. 



^base 


dtso. 1 

L. — ^ s2 

dC^ 


where. 



CHAPTER 3 


EVALUATION OF SUBLANCE MODELS FOR END POINT 
CARBON AND TEMPERATURE 


3.0 Introduction 

Various models for end point control of carbon and temperature have already been 
discussed in Chapter 2. Operational data for 300 tons and 210 ton converters used in the 
present work are described in section 3.1. Results obtained for some of these models are 
presented in section 3 . 2 . So as to verify their applicability to different converters and then 
find out the best model for a given situation. Comparative performance of different models is 
discussed in section 3.3. 


3.1 Operational Data 

The data collected from 4 plants is analyzed in this work . The first set 
comprises of 301 heats made in a 300 ton converter using an oxygen flow rate of 600 m /min. 
In plant data set 1 (see Appendix 1) the variables included are ; hot metal weight (tons), scrap 
(tons), hot metal silicon (milli % ), oxygen blown in second blow (O 22 ), temperature at 
sublance measurement (Tsl in Kelvin), carbon at sublance measurement (Csl in milli 
percent), end point temperature (Tf in Kelvin), end point carbon (Cf in milli percent), lime 
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added in first and second blow (kg), return slag in first and second blow (kg), dolomite 
added in first and second blow (kg), ore added as a coolant (kg), dolomite special added (kg), 
lance height (m) and lance life (m). 

The plant 2 data of two different converters (but of same capacity) comprises of 1744 
heats. The variables included in plant 2 data set (see Appendix 2) are hot metal weight (kg), 
silicon in hot metal (in milli % ), lime added in first blow (kg), carbon and sublance 
measurement (Csl in milli % ), end point carbon (Cf in milli % ), temperature at sublance 
measurement (Tsl in Kelvin ), end point temperature (Tf in Kelvin), oxygen blown in second 
blow (m^), ore added as a coolant (kg). The plant data set 3 (see Appendix 3 ), comprises of 
748 heats of a partial campaign of 3 15 ton converter in which oxygen blowing rate was 900 
mVmin. The variables included in set 3 are hot metal weight (kg), silicon in hot metal (in 
milli % ), lime added in first blow (kg), carbon at sublance measurement (Csl in milli % ), 
end point carbon (Cr in milli % ), temperature at sublance measurement (Tsl in Kelvin ), end 
point temperature (Tf in Kelvin), oxygen blown in second blow (m^), ore added as a coolant 
(kg). 


The fourth set, plant 4 data set, (see Appendix 4), comprises of 2235 heats of a 210 
ton converter, in which oxygen flow rate is 600 m^/min.. The variables included in plant 4 
data are : hot metal weight (kg), silicon in hot metal (in milli % ), lime added in first blow 
(kg), carbon at sublance measurement (Csl in milli % ), end point carbon (Cf in milli % ), 
temperature at sublance measurement (Tsl in Kelvin ), end point temperature (Tf in Kelvin), 
oxygen blown in second blow (m^), ore added as a coolant (kg). 
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Some special features of different sets of data are as follows. Set 1 includes 
additional data on lance life and lance height. Set 2 comprises of a large number of heats thus 
making it possible to test the models over a full campaign. Set 3 is for a converter in which 
very high flow rate of oxygen (900 mVmin ) is used and Set 4 is for a converter in which 
bottom stirring was absent. 

3.2 Results 

3 . 2.1 Application of the Model by Chou et al. 

According to Chou et al.[2], the decarburization rate equation as given in chapter 2, (Eq 
(2.24)) is 

4%C] -(2X1 

dt Ulx\(r\\-Q5Xco)W, 

A sample calculation can be done for a converter with 
Oxygen flow rate = Qo = 900 Nm^/min 

Weight of bath = Wt =300tons = 300x lO^gm 

Relative component of CO in (CO+CO 2 ) = Xco = 0.9947 at [%C] = 0.3 
On integrating Eq (2.24) in the limit 4 % C to 0.3 %C, we get 

_ Qo (1 ^Inr^ ^ 

^ ° l.S7xl0-\l-0.5XcoWi' 


(3.1) 
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Assume the gas used is pure i.e., xinr = 0, then 


(0.3-4) = 


-900(1) ^ 

I.87a: 10'-'(1 -0.5x0.9947)x300xl0® ’ 


Time required to reduced carbon content from 4% to 0.3% becomes 


t = 11.59 min. 

Equation(2.27) can be written as 

-4%C] ^ ro/,ci 
dt V 


On integrating the above equation in the limit of 0.3 to 0.04 %C, we get 


0.04 = 0.3 exp (-kd. A/V. (t- 1 1.593)) 


(3.2) 


Time required for reducing carbon content from 0.3% to 0.04% can now be calculated from 
the above Eq (3.2), provided kd.(AA'') is known. 

From Eq (2.24) at [%C] = 0.3 


£/[%C] -900 

dt ~ 1.87;cl0"*(l - 0.5x0.9947)x300xl0® 


c/[%C] 

dt 


-0.3192 
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Substituting in Eq. (2.27) 


-0:2128 = — ^[0.31 
V 


Or kd.(A/V) = 1.064 


Substituting this value in Eq (3.2) 


4%C] 

dt 


-1.064 [%C] 


On integrating the above equation in the limit of 0.3 to 0.04 %C, we get 

[%C] = 0.3 exp (-1.064 (t - 1 1.593)) 

Or (t- 11.593)= 1.89 min = 113.62 sec 

The total blowing time for second and third stages of blow (25% to 100%) is thus 

t= 11.593 + 1.89 
Or t= 13.48 min 

In the first stage of blow (approximately 3 minutes ), primarily silicon is oxidized. 
Total blowing time is thus 13.48+3 =16 min. 

This is approximately same as observed in an actual converter. It proves that the 
model developed by Ghou et al. can be tuned to specific operating conditions in a given 
converter to predict end point carbon content in an actual converter. The model developed in 
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present work (as discussed in Section 3.3) is actually an improvement over the model of 
Chou et al. 

3.2.2 Application of Exponential Model 

As explained in Chapter (2), the decarburization rate is given by Eq (2.29) as 




where Xp’ is dependent upon operational variables. 

After several preliminary evaluations, the operational variables chosen for the 
regression in this work are 


Xi = 

hot ratio (wt. of bath/wt. of scrap) 

Xz = 

bath temperature by sensor lance (°K) 

Xs = 

slag volume (m^) 

X4 = 

second blow addition (kg) 

Xs = 

ore added during end blow as coolant (kg) 

Xfi = 

lance life (m) 

Xt = 

lance height (m ) 


In Model No.(l) operational variables Xi to X? are considered. In Model No. (2) 
operational variables Xi to Xe are considered. Lance height (X?) is used in sublance model 
described later. 


Model No. I : 
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By using stepwise regression, which selects only the significant variables significant in 
affecting the dependent variable 'Cp’, only Xi and X 5 were selected in the regression. It 
is, therefore, assumed that for the operational data chosen only these two variables are 
affecting the dependent variable (Cp). 

The value of Cp can be calculated from the Eq. (2.30) and the end point carbon can be 
calculated from Eq. (2.3 1). A sample calculation for one heat is shown below. 

Sample Calculntion 


Step I ; From the plant data 'Cp ‘ is calculated by solving the non-linear Eq. (2.30) 
by iterative process. 


Aa 


(.c.-cp) 




■exp 


In- 


C -C . 
C. -C„ 


C..-C 


•exp^ 


SL 


C -C 

^'oa ^ fSL 

c 


Cp-c; 


here ^ = 0.214 for AO 2 = 4.877172, Co =10 Cta = 61, Csl = 257 
Cp value is found out to be = 0.6069 


Step 2 : With this ' Cp‘ and the operational variables calculated from the plant data, 

regression is done according to the Eq. (2.32) 




i=\ 


For these two selected variables, the regression coefficients are found to be 


Xi = hot ratio = 6.0 ; ai 


0.0187208 
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Xj = ore added = 1004 ; aj = -0.0000532 

Regression constant = ao = 0.484089 

Step 3 : Cp is calculated from Eq. (2.32) by substituting regression coefficients 

i.e. Cp =(0.484049) + (0.0 187208)(6.0)+(-0.0000532)(1004) 

= 0.543001 


Step 4 : Now, Cp is substituted in Eq. (2.3 1) to estimate the end point carbon. 


Q. = (Q-r)- 


R 


C - C, 


l_exp~2 — ^ + i?.exp- 


C„ 


C,-C 




where. 


R = exp 


Cp-c; J 


(2.35) 


For 

Co = 10 

Cp = 0.543001 

CsL = 257 

A 0’2 = 4.877172 

Cfa is found out to be 54.2 milli percent. Actual end point carbon content for this heat 
is 61 miili percent. 

As explained in the sample calculation, the end point carbon is estimated for 
all the heats. The plot of actual end point carbon values versus estimated end point carbon 
values (from the plant data) is shown in Fig. (2) . 
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The results of this regression for actual versus estimated end point carbon values are 
as follows (also summarized in Table ( 1 ) ). 


For the regression curve 

Befor e ± 2 O’ 

After ± 2 0 

y = mx 

II 

U.) 

0 

n = 284 


m = 0.9833 

m= 1.00528 


0 = 14.002 

a = 11.4872 


R = 0.6612 

R = 0.73016 


F = 6885 

F = 9404 

y = mx+C 

n = 301 

n = 285 


m = 0.5981 

m = 0.5958 


C = 27.3968 

C = 27.002 


0 = 12.041 

0=10.1623 


R = 0.6612 

R = 0.71 


F = 232 

F = 287 


In Model No. 1 the variable slag volume was replaced with basicity to test whether 
basicity has an effect on Cp or not. But again the same variables Xi (hot ratio) and X 5 (ore 
added during end blow as coolant) are selected. 


Model No. 2: 


Lance life is taken as additional variable to estimate Cp. Dependent and independent 
variables tried are as follows; 


Dependent 
variable in the 
regression 

Independent variables considered 
for the regression 

Selected variables by stepwise 

regression and their corresponding 
coefficients 

Cp 

Xi = Lance life 

X2 = Temperature at sublance 

X3 = Second blow additions 

X4 = Slag volume 

Xj = Ore added during end blow 

X 6 = Hot ratio 

Xi =lance life (ai = -0.0002785) 

X 5 = ore added (as = -0.0000572) 

X 6 = Hot ratio (a^ = 0.01855860) 
Regression constant (ao = 0.509299) 






Estimated end point carbon C W IL' 7. ) Estimated end point carbon C W Ua' 


f(x)=0.5981*x+27.3968 


a=12.0415 

R=0.6612 

n=301 






0,0 o 


^ A O o'^ 

o^ % o 

Oo^ #oo o 

5^ o.9’o ' ^ 

Ivlw „» « 


:o o c 


t 


Actual end point carbon 

f 'TA^iA. 


After removing + 2 ct points 


f(x)=0.5958*x+27.002 


a=l 0.1623 


R=0.71 


oi> ^ % 


100 - n=285 o 0 0^0 

Oo ^oo ^ 

40-"o\%<^ 

o. <> ^ I I 


Actual end point carbon 


^ov\>u1Aa /, ^ 


Fig. 2 Actual end point carbon versus Estimated end point carbon 
for Model No. 
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With the above selected variables, Cp is calculated by substituting the coefficients in 
Eq (2.32), as follows. 

Cp=0.509299-0.0002785( 30 )-0.0000572 (1004 )+0.01855860( 6.0 ) 

=0.55487 

The Cp values are substituted in Eq (2.35) to estimate end point carbon as explained 
earlier. The results of regression are as follows (also summarized in Table No. ( 1 ) ) 


For the regression curve 

Before ±2 a 

After ± 2 a 

Y = mx 

n = 301 

n = 292 


m = 0.98808 

m = 0.9913 


a = 13.7495 

CT= 12.1226 


R = 0.66998 

R = 0.71026 


F = 7102 

F = 8724 

Y = mx+C 

n = 301 

n = 287 


m = 0.5782 

m = 0.5809 


C = 27.3066 

C = 28.0326 


a = 11.7602 

a = 9.9284 


R = 0.67 

R = 0.7081 


F = 243 

F = 286 


The standard deviation (a ) in estimating end point carbon is 9.9284 milli percent. 
This estimation is better than the earlier one (Model No. 1 ) i.e. without the lance life for 
which a obtained is 10.1623. 


Model No. 3 ; 


An attempt has been made to reduce the standard deviation in estimating the end point 


carbon by incorporating mean values of operational variables in MLR equation(2.32). Thus 
the MLR equation is modified as 
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C,=a + -£a,iX,-X,) (3.3) 

1 = 1 

Xj =mean value of operational variable X, 

According to the above equation, Cp is calculated and then substituted in Eq (2.35) to 
estimate the end point carbon. 


Dependent variable in the 
regression 

Selected variables by stepwise regression and their 
corresponding coefficients 

Cp 

Xi =ore added-mean of ore added (ai = -0.0000532) 

X 2 = (hot ratio-mean of hotratio) (a .2 = 0.01 87208) 
Regression constant (ao = 0.5556) 


With these selected variables Cp value is calculated by substituting the coefficients in 
MLR Eq (3.3) as follows. 

Cp=0.5556-0.0000532(409.475)+0.0 187208(0.489400) 

=0.54297 

These Cp values are substituted in Eq (2.35) to estimate the end point carbon values. 
The regression between estimated end point carbon and actual end point carbon has been 
done. The plot of actual end point carbon versus estimated end point carbon has been shown 
in Fig (3). The details are as follows. 





Estimated end point carbon ^ 7-) Estimated end point carbon 



40 60 80 100 120 140 

Actual end point carbon 


After removing ± 2'6’points 



40 60 80 100 120 140 

Actual end point c^bon 

(_^w\^ IAa 7 ^ 


Fig. 3. Actual end point carbon versusEstimated end point carbon for Model No. (3) 
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For the regression curve 

Before ± 2 a 

After ± 2 a 

Y = mx 

n = 301 

n = 286 


m = 0.96853 

m = 0.95147 


a = 13.8538 

G= 11.2682 


R = 0.661 19 

R = 0.7349 


F = 6865 

F = 9373 

Y = mx+C 

n = 301 

n = 284 


m = 0.7309 

m = 0.6916 


C = 16.6945 

C = 17.6208 


a= 13.31 

a = 10.3263 


R = 0.6612 

R = 0.7301 


F = 232 

F = 321 


As shown above, standard deviation is decreased from 11.4872 (In Model No'2 
for y=mx) to 11.2682 in Model No. 3 by incorporating mean values of operational variables 
in MLR equation. It is therefore proven that by incorporating mean values of operational 
variables in the regression a slightly better result can be obtained. 


3.2.3 Sublance Model Based on Mass transport Control 


On integrating Eq (2.5) from Csl to Ct , we get 


-In 


r c ^ 

\^SLj 


, A 
= k^.—l 
^ V 


Or 


•In 


\^SLy/ 


K.t 


where 


K=k'- 


(3.4) 


(3.5) 


Ct =end point carbon content 
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CsL =carbon content at sublance measurement 
and t =time to reduce from C, to Csi 

The term kc-A/V is actually called ‘the capacity coefficient of mass transfer ‘ and is 
usually represented by K. This K is specific to a converter practice at a particular plant. 

The value of K may depend upon lance nozzle design, lance height, slag volume and 
composition, coolant additions in the converter, etc. In the present work, owing to the 
complicated nature of process we will be concerned with the overall value of K and not the 
individual values of the terms kc, .A and Y. This is because the phenomena of recirculation 
of metal in jet impact zone, droplet disintegration (due to the evolution of CO gas ) as well as 
droplet coalescence during their passage in slag is still not well understood. Droplet 
disintegrates into hundreds of minute droplets (sparking ), as visually observed during entire 
converter operation. Coalescence of droplets would also depend upon the inter-collision of 
droplets during their zig zag passage through the foamy slag. 

In the present work two methods have been adopted to incorporate the effect of 
operational conditions on decarburization. 

Method 1 : 


It is assumed that ‘K’ is directly dependent upon the different operational variables. 
Multiple linear regression with different operational variables is first done as follows. 

First ‘K’ is estimated for each heat from the equation 
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K = 


-In 


'SL'' 



The value ot 'K' tor each heat is then related to the operational data by 


j=i 


where. 


Xi = operational variables 
ai = regression coefficients 
ao = regression constant 

The calculated value of ‘K‘ is substituted back in 


C, =Qi.exp(-A:./) 


to estimate the end point carbon. 


Method 2 : 


In this method , Eq (3.12) is modified as 



1 


C n ^ 

In 

= K.t + 



J 


1 -=! y 
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and the regression is done for all the heats to obtain values of K, ao and aj. Then Ct is 
estimated from 


C, = C5^.exp[-^./ 


a,+^a..x, ] 
V .=! J 


.-■> . I‘ ? ' 


In both the above mentioned methods, selection of operational variables is a critical step. For 
this “stepwise regression iTiethod“ is used . This technique selects only the variables which 
are significant in affecting the dependent variable. 

In the above equations the variable ‘ time ‘ for each heat is calculated in two 

ways. 


time = 


oxygen. blown m.the. sec and. blow 
oxygen. flow.rate 


In an alternative approach the time is normalized for total iron content of bath by assuming 
that scrap contains approximately 100% Fe and hot metal contains approximately 92% Fe. 

^ oxy gen. blown. in Ihe. second. blow. X. Total Wt. of .bath 

oxy gen. flow. rate. X t. of .bath. X 0.92 + Wt. of scrap) 

These two methods are applied to plant (1) data set first. Then after deciding the method 
which gives a better result, that method is applied to the Data sets 2-4. 
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3.2.3.1 Application of Method I for plant Dataset 1 : 


As summarized below, K is related to six operational variables. The selected 
variables are only basicity (Xi) and ore added during end blow (X5), as shown below. 


Model No. 4 : 


Dependent 
variable in the 
rearession 

Independent variables considered 
for the regression 

Selected variables by stepwise regression 
and their corresponding coefficients 

K 

Xi = Temperature at sublance 

X 2 = Second blow additions 

X3 = Basicity 

X 4 = Ore added during end blow 

Xft = Hot ratio 

Xi =Basicity (ai = -0.000215759) 

X5 = ore added (as = -0.000000635) 
Regression constant (ao = 0. 10422) 


With these selected variables ‘ K ‘ values are calculated by substituting the 
coefficients listed above as follows. 

l=z\ 

K=0.10422-0.000215759(3.9990)-K).000000635(1004) 

=0.10399 


These’ K ‘ values are substituted in Eq (3.8) to estimate the end point carbon values. 
The regression between estimated end point carbon and actual end point carbon has been 
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done. The results are as follows. The plot of actual end point carbon versus estimated end 
point carbon has been shown in Fig (4). 


For the regression curve 

before ±2 a 

after ± 2 a 

Y = mx 

n = 301 
m = 0.496 
a= 13.33 

R = 0.6627 

F = 7046 

N = 288 
m = 0.96759 

a= 11.43 

R = 0.73 15 

F = 9130 


The result obtained here i.e. by sublance model ( a =13.33 mlli % ) is better than 
obtained by exponential model ( a =14.002 milli % in Model No.l) for the same data set. It 
may be noted that only 3 data points are discarded for data points lying outside ± 2 a limits. 


Model No . 5 : 


If lance life is included in the regression, it is found that a slight improvement in 


results can be obtained. 


Dependent 
variable in the 
regression 

Independent variables considered for 
the regression 

Selected variables by stepwise 
regression and their corresponding 
coefficients 

K 

Xi=lance life 

X2 = Temperature at sublance 

X3 = Second blow additions 

X4 = Basicity 

X5 = Ore added during end blow 

Xe = Hot ratio 

Xi =lance life (ai = 0.000003467) 

Xs = ore added (as = -0.000000683) 
X 4 =basicity (a4==-0. 000219421) 
Regression constant (ao = 0.10135) 










Bsciilialcd end poiiii carbtm AA^ / J lisiinialcd end point carbon 
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With these selected variables ‘ K ' values are calculated by substituting the coefficients in 
MLR Eq (3.6) as follows. 

1=1 

K=0. 1 0 1 3 5-0.000003467(30)+0.000000683 ( 1 004)-0.00021 941 2(3.9990) 

=0.10105 

These’ K ‘ values are substituted in Eq (3.8) to estimate the end point carbon values. 
The regression between estimated end point carbon and actual end point carbon has been 
done. The details are as follows. 


For the regression curve 

Before ± 2 a 

Y = mx 

n = 301 
m = 0.95390 
a= 13.173 

R = 0.6688 

F = 7211 


Model No. 6 : 


If the variable basicity is replaced with slag volume to see the affect the result 


obtained is as follows. 
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Dependent 
variable in 
the 

regression 

Independent variables 

considered for the regression 

Selected variables by stepwise 
regression and their corresponding 
coefficients 

Regression 
of estimated 
.Vs. Actual 
end point 

carbon. (y=m 

X) 

K 

Xi=lance life 

X 2 = Temperature at sublance 

X3 = Second blow additions 

X4 = slag volume 

X 5 = Ore added during end 
blow 

Xe = Hot ratio 

X, =lance life (ai = 0.009085) 

X 5 = ore added (as = 0.000003288) 
Regression constant (ao = 0.009085) 

m=0.95329 
a =13.522 
R=0.6497 
F=6834 


As shown above , ‘ a ‘ is increased from 13.173 to 13.522 (when compared with 
Model No. (5 )). Thus incorporation of ‘ basicity ‘ as independent variable gives better 
result than that of’ slag volume’. In subsequent Models only basicity will be used to estimate 
end point carbon but not the slag volume. 

3.2.3.2 Annlicntion of Method 2 for plant Data set 1 : 

In this method, the term ‘-ln( Q/ Csl)‘ is related to the following operational 

variables. 


Model No. 7 : 


Dependent 
variable in the 
regression 

Independent variables considered 
for the regression 

Selected variables by stepwise regression 
and their corresponding coefficients 

-ln( Ct/ Csl) 

Xi=time 

X 2 = Temperature at sublance 

Xj = Second blow additions 

X 4 = Basicity 

Xs = Ore added during end blow 

Xe = Hot ratio 

X, =time (K= 0.007827967) 

Xs = ore added (as = 0.000 1143 56) 
X 4 =basicity (a4=-0.04003613 1) 

Regression constant (a<, = 0.473776) 
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With these selected variables * - In 




is calculated by substituting the 


coefficients K, ao and ai in MLR Eq (3.9) as follows. 


-In 






= KJ + 


( ^ 

V i=l ^ 


-In 


K^sl j 


=0.007827967(152)+(0.473776+0.000114356(1004)-0.040036131(3.99)) 
= 1.618 


From these 


■In 


A 

C 

y^sLj 


values, ‘ Ct ‘ values are calculated as follows. 


-In 





C, 


SL J 


= 1.618 


From Eq (3.10) 


Ct=CsL exp (-1.618) 

= 257 .exp (-1.618) 

=50.926 

The plot of actual end point carbon versus estimated end point carbon has been 
shown in Fig (5).The estimated end point carbon values are then regressed with actual ones 


and the results of regression are as follows. 


For the regression curve 

Before ± 2 ct 

After ± 2 CT 

Y = mx 

n = 301 

n = 289 


m = 0.95311 

m = 0.97334 


CT = 12.9347 

CT= 11.027 


R = 0.66582 

R = 0.73 165 


F = 7466 

F = 9867 

Y = mx+C 

o 

II 

d 

n = 282 


m = 0.496 

m = 0.497 


C = 32.0134 

C = 31.4141 


CT = 9.8615 

CT= 8.3363 


R = 0.6658 

R = 0.7138 


F = 238 

F = 290 
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It can be observed that the ‘o ‘ obtained in this model (12.9347) is less than, the ‘ a 
‘obtained (13.33) by the tirst method (Model No. 4 ). It can therefore be concluded that 
method ( 2 ) gives a better result than method ( 1 ). 


Model No. 8 : 


If lance height is included in the regression, we get a slightly better prediction. 


Dependent 
variable in the 
regression 

Independent variables considered 
for the regression 

Selected variables by stepwise regression 
and their corresponding coefficients 

-ln( C,/ Csl) 

Xi=tiine 

X 2 = Temperature at sublance 

X 3 = Second blow additions 
.X 4 = Basicity 

X 5 = Ore added during end blow 

X(, = Hot ratio 

X 7 =lanceheight 

Xi=time (K= 0.0078753) 

X 5 = ore added (as = 0 . 0001102 ) 
X 4 =basicity (a 4 =- 0 . 0466949) 

X 7 =lanceheight (a7=-0.0066976) 
Regression constant (ao = 2.013246) 


With these selected variables '-In 





is calculated by substituting the 


coefficients K, ao and ai in IVILR Eq (3.9) as follows. 



C, ] 

[CsJ 


f 

In 

= KJ + j 




1 1=1 J 


-In 




0.0078753(152)+(2.013246-0.0066976(227)+0.0001 102(1004)- 


= 1.61426 


0.04669496(3.99)) 





a =9.8615 
R=0.6658 
n=301 


f(x)=0.4964=*‘x+32.03 14 


o<» 


o 

o. 


4o%Oo^ oO O 0 O 


60 80 100 120 
Actual end point carbon 

CWU; 7.) 


After removing ± 2 O’- points 


^ O =8.3363 

f R=0.7138 

O 120 - n=282 


f(x)=0.497*x^3 1.4141 


o o Oo . . - o 'o ^ 
o$ . ..O’ "' 


<» 








o- 


40 


60 80 100 120 
Actual end point carbon 

C'YvuHCV.) 


Fig 5. Actual end point carbon versusEstimated end point carbon for Model No. (7) 
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From these 






values. 


Ct” values are calculated as follows. 


-In 


'll' 

SLJ 


1.61426 


From Eq (3.10) 


C, =CsL.exp (-1.61426) 
= 257. exp(- 1.6 1426) 
= 51.15252 


The estimated end point carbon values are regressed with actual ones. The results of 
regression are as follows. 


For the regression cuive 

Before ±2 a 

After ±2 a 

Y = mx 

n = 301 

n = 261 


m = 0.95856 

m = 0.98009 


a= 11.816 

a= 10.012 


R = 0.73992 

R = 0.79227 


F = 8226 

F = 11022 

Y = mx+C 

n = 272 

n = 259 


m = 0.5692 

m = 0.5977 


C = 27.3729 

C = 25.5638 


a = 9.333 

a = 8.092 


R = 0.7399 

R = 0.7784 


F = 326 

F = 395 
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It can be observed that the a ‘ obtained (8.092) with lance height as a variable is 
less than, the ‘ a ‘ obtained (8.3363) in Model No. (7) i.e. without lance height. So, it is 
concluded that with the known lance height values, better estimation can be done. 

Model No. 9 : 


As a further step to improve the estimation of end point carbon, the equilibrium 
concentration of carbon (C.-) is incorporated in the LHS of the MLR equation i.e. the MLR 
Eq(3.9) is modified as 


f c -C 'l 


f " 


= K.t + 

ao+^a-.x,. 

^ <=i ; 

c - c 

'-c-J 


and the regression is done for ail the heats to obtain values of K, ao and a,-. Then ‘Q’ is 
estimated from 


Q = C, + < 


(C^-C,)exp 


K.t + 






1=1 


i3,!4i 


The term - In 


yC^-C.j 


is related to the following operational variables. 
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Dependent 
variable in the 
regression 

Independent variables considered for 
the regression 

Selected variables by stepwise 
regression and their corresponding 
coefficients 

( C - C 
-In — ^ ^ 

Iq-cJ 

Xi=time 

X 2 = Temperature at sublance 

X3 = Second blow additions 

X 4 = Basicity 

X5 = Ore added during end blow 

X* = Hot ratio 

X, =time (K= 0.00855417) 

X5 = ore added (a.. = 0.000 1 18546) 
X 4 =basicity (a 4 =- 0 . 04342846) 

Regression constant (a<j = 0.505659) 


With these selected variables '-In 


^ C,-C, ^ 


V ^ SL 


■c,. 


• is calculated by substituting the 


coefficients K, ao and ai in MLR Eq (3. 13) as follows. 


-In 


( C -C ^ 


( 3 ^ 

'“f c- 

— K .t + 




V -=i y 


=0.00855417(15 



1.75163 


From Eq (3.14) 




= 10 +[(257-10) .exp(-1.75163)] 

=52.8521 

The estimated end point carbon values are regressed with actual ones. The plots are 


shown in Fig (6). The results of regression are as follows. 
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For the regression curve 

Before ± 2 o 

After ±2 0 ' 

Y = mx 

n = 301 

n = 286 


m = 0.94268 

m = 0.96702 


CT= 12.6652 

o= 10.5013 


R = 0.66826 

R = 0.73212 


F = 7618 

F=10511 

Y = mx+C 

II 

O 

n = 286 


m = 0.4428 

m =0.45 


C = 35.0092 

C = 34.1851 


CT = 8.7458 

0 = 7.613 


R = 0.6683 

R = 0.712 


II 

u. 

F = 291 


It can be observed that the ' a ' obtained here (7.613) is less than, the ‘ cr ‘obtained (8.3363) 
in Model No. (7) for the same set of data. Thus, it is proven that better estimation of end point 


carbon can be done by using the term -In 


( c,-cA 


i.e. considering the effect of Ce 


instead of using the term - In 




Model No. 10 : 


If lance height is included in the regression, we obtain a better result (as explained in 


Model No. 8 ) 


For the regression curve 

Befor e ±2 o 

After ± 2 o 

Y = mx 

n = 272 

n = 263 


m = 0.94839 

m = 0.96813 


0=11.61 

o = 10.0727 


R = 0.74384 

R = 0.78413 


F = 8341 

F= 10739 

Y = mx+C 

n = 272 

n = 262 


m = 0.517 

m = 0.5284 


C = 30.3324 

C = 29.3979 


0 = 8.3772 

o = 7.4204 


R = 0.7438 

R = 0.7728 


F = 334 

F = 385 





m 


a -8.7458 

R-0.6683 

n=301 

f(x)=0.4428*x-H35.0092 





* 
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Fig. 6. Actual end point carbon versus Estimated end point carbon for Model No. (9) 
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3.2.3.3 Results for Data sets of plant (2) , plant (3) and plant(4) : 

It has been shown earlier that Model No.(9) based on Eq.(3.13), (modified equation 
suggested in method (2)) gives the best result. This method is now applied to estimate the 
end point carbon. 

Results obtained for plant data sets (2) and (3) are summarized in Table No.(l) 
under Model No. (II), Model No.(12), respectively. It may be noted that the lance height 
details are not available in data sets (2) and (3). Otherwise, by using lance height variable in 
the regression, a lower standard deviation in estimating the end point carbon could have been 
achieved. It can also be observed that ‘ a ‘ for plant data set (2) is lower than for plant data 
set (1). This is because, the plant ( 2 ) data set has more number of data points(1744 data Vs. 
301 data). 

The results of carbon prediction for plant 4 data set are summarized in Model 
No. (13) and Model No. (14), in Table (1). In this set of plant data, carbon is predicted for 
only for those heats in which end point carbon is less than 0.06 %C. The standard deviation 
for plant 4 data set is highest because of the fact that bottom stirring was absent. 

3.2.3.4 Application of Temperature Model 


As explained in Chapter (2), N.Bessho et al. have given a model, Eq(2.33), to 
estimate end point temperature 


t , = t ^+ f ( a . o,)+q 


yCj- 


c 


+R 


SL J 
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Fig 7. Actual end point carbon versus Estimated end point carbon for Model No. 1 1 
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This model is applied in the present work to estimate the end point temperature by taking 
different combinations of operational variables which may effect the end point temperature. 
Various alternatives have been tried to estimate the end point temperature with the minimum 
possible standard deviation. The results of all the Models for plant data set 1 are summarized 
in Model Nos. (15),(16),(17), in Table (1). 

In the above Eq(2.33), the end point carbon Cf is not known a priori at the 
time of sublance measurement. So, it is not correct to use this information to predict the end 
point temperature though it gives a good result. Therefore the equation is modified as 


= 7;,,+P(A.O,)+fJ 


^ 1 ^ 




+R 


(3.15) 


This equation was applied to plant data set 1 to predict end point temperature and the results 
are summarized in Model No.(18),(19),(20),(21),(22), in Table (1). 

The Eq(2.33) was applied to predict end point temperature by substituting the 
predicted end point carbon (Cr) from the sublance carbon model in Model No. (23). The 
results are summarized in Table (1). 

Similarly for plant data sets (2) and (3), Eq(3.15) is applied to predict the end point 
temperature. The results are summarized in Model No. (24) and Model No.(25), respectively, 


in Table (1). 

For plant data set (4) , the end point temperature was also estimated from the 
equation (2.48) which was given by Byun et al.(as explained in Chapter (2)). 




i. 

C. 


'SlJ 


+5, .AO, 
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The results are summarized in Model No.(26) in Table (1). 

The earlier Model given by Bessho et al., Eq(2.33) was used to estimate the end point 
temperature by substituting the end point carbon estimated from the sublance carbon Model 
No.(13). The results are summarized in Model No. (27) in Table (1). 


3.3 Discussion 


Fundamental aspects of actual mechanism of decarburization (i.e. droplet 
versus bulk decarburization) are discussed in section 3.3.1. This is followed by comparative 
evaluation of the results obtained by using models for different sets of plant data in section 

o O 

3.3.1 Estimation of mass transfer coefficient and droplet decarburization 
versus bulk metal decarburization 

The two possible mechanisms of decarburization are 

• Decarburization via droplets ejected from the jet impact zone and 

• Decarburization at the interface of bulk metal and slag and in the recirculation zone 
below the jet. 
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Fig 10. Actual end point temperature versus estimated end point temperature for 


Model No. 15 
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Fig 12. Actual end point temperature versus estimated end point temperature for 


Model No.24 
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The metal droplets are produced due to the jet impact and the shearing action of the 
gas flow from the impact region when the jet strikes the metal surface and the product gases 
are deflected upwards. 

The limiting diameter of the droplets (i.e. largest diameter droplets, diimit) is a 
function of the pressure at which oxygen is blown, lance height, lance inclination and nozzle 
diameter. Consider a 300 ton converter for which 
Flow rate of oxygen (F ()2 ) = 915 m Vmin, STP 

Stagnant Pressure P„ = 10.2 atm 

Throat dia dt = 46,6 mm 

Exit pressure P,- = 1.0 atm 

Exit dia = 58.0 mm 

Nozzle angle 8t = 12“ 

No. of notszies n = 5 

Tbuik = 1973 K 

Tho,sp.« = 2200 K 

Lance height == x = 2 m 

The expression to calculate diinm is given by [12] 

= 5.513xl0“'[10'’(— )'-PJI.27A-l}Coi-^]‘"'^OTOT 
X 

= 5.5l3rf0->[10'-(25f5i)=xlxjl.27(^)-nCav(12)]'”‘ (3.16) 

X 1 

1130.95 

At a lance height of 2 m (x = 2) , diimii = 212.5 mm (or = 21.25 cm) 

Average diameter of droplets is given by [12] (davg) “ 0-2 diimit - ^2) 


= 0.2 (21.25) = 4.25 cm 
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Thus average radius of a droplet = 4.25/2 = 2. 125 cm 

If diffusion coefficient of carbon is 3.24x10^ cmVsec then from Eq. (2.6) the mass transfer 
coefficient kc can be calculated as 

kc = 5D/a = 5 X 3.24x10'^) / (2. 125) 

= 7.623x10'’ m/sec 


From Eq (2.5) the decarburization rate is 


-<j[n 

dt 




On integrating the above equation, we obtain 


r 

■ln( ^ ' 


C -C 

^ SL ^ e 


C. , , A , , 

'-) = k,. —t( — ) 

f / C T r ^ XT 


V Q-q 


(3.18) 


where A = surface area 
V = volume 


If the equilibrium carbon concentration 'Ce’ is assumed to be very small then the 


equation (3.18) becomes 






SL 


A 

V~ {A [2)70-^ 


2'=_L_ = 1.4117 cm‘‘ 
r 2.125 

= 141.17 m'^ 


(3.19) 


where 
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From Eq (3.19) the time required to decrease the carbon concentration from 0.3% to 
0.04% can be calculated as 

-In (0.04/ 0.3) = 7.623 X 10'^ x 141.17 xt 
Or t = 187 sec. 

According to the procedure suggested in [12], iron conversion (i.e., total mass of iron 
droplets ejected from impact zone (tons/min) ) can be calculated as follows. 

The dynamic impact pressure (Pa..,;) at the jet impact is given by an empirical 
relation[12]. 


= 23 0(—)'"" (3.20) 

Po dt 


where, 

Pd.x= Po-Px 

X = distance from nozzle exit 

The functional relationship between iron conversion and Weber number (We) is [12] 

Mg = f(We) (3-21) 

where the Weber number'is defined as 

We = _ "^Pdx 


( 3 . 22 ) 
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where, 

p, = density of liquid bath (7000 kg/m’^ ) 

Pg = density of gas 

cr, = surface tension of liquid ( 1 .7 N/m) 

Pd,x = Po-Px 

Uojj = center line velocity of gas at distance x from nozzle tip 

In order to apply Eq{3.6) to a real system, the droplet generation per unit volume of 
blown gas must be considered. The plot of Mg/Fo 2 versus Weber number is shown in Fig. 
(16). Using this graph, the best fit polynomial has been obtained as follows; 


Mg/Fo 2 = 9,H02 - 0.8816 x (We) + 0.02368 (We)^ - 0.1415 x 10'^ (We)^ + 0.1572 x lO^^ 
(We)"^ +0.5402 x 1 0'“’ (We)'' (3.23) 

On substituting Eq (3.20) in Eq (3.22), we get 

(3.24) 

1.7 N/m 

Substituting this value in Eq. (3.23) 

Mg/Fo2 = 56.65 kg/m^ 

Or Mg = 56.65 x 915 mVmin = 51832 kg/min 


y, - 2.0 x230.x-/> ^ .r 
ipiga,) dl 

since x =2 m, g = 9.8 1 m/s^ ,Po =10.2 atm, dt =0.0466 m, pi=7000 kg/m^ , ai = 
we get We =167.88 



Mg/F02 
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S8 


= 51.832 toii/min (or s: 50 tons/min) 

Thus, the mass of droplets continuously circulated through slag during 
decarburization period is approximately 50 tons/min and only this metal would be subjected 
to decarburization via droplets (i.e. approximately one sixth of total metal). Additional 
decarburization can be conceived to occur either by assuming an enhanced mass transfer 
coefficient (by a factor ot 6) within the droplet and/or decarburization at the interface of bulk 
metal and slag, and/or decarburization in the impact zone due to high recirculation rate of 
metal in the zone. 

As far as bulk metal decarburization is concerned, the nominal area of metal bath for a 
300 ton converter is approximatelv 30 m". Thus 


r 300x1000 
7000 

The value of (A/V) for droplet is 141.17 in'* where as (A/V) for bulk metal is 0.7, which is 
lower approximately by a factor of 200. It shows that the contribution of decarburization at 
the interface of bulk metal and slag is expected to be very small. 

The diameter of jet impact zone is approximately 1 m. For a lance tip with 5 
nozzles, the total area of impact will be 5x tc x (0.5)^ =3.92 m^ . Suppose that entire volume 
of metal is recirculated then V=42.88 m^ . Now for Model No. 4 the value of capacity 
coefficient of mass transfer is approximately 0.1 s ' . Hence mass transfer coefficient can be 


computed as 
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Or k,^K}L = 0.\x^mls 
A 3.92 

= 1.09 m/s 

Since nitiss tuiist'er coefficient is directly proportional to fluid velocity, the 
average recirculation velocity of metal is of the order of 1 m/s. Such a magnitude has been 
predicted from fluid flow computations in the impact zone[12] . It therefore indicates that 
most of decarburization can occur in impact zone. 


3.3.2 Comparative evaluation of models for carbon and temperature 

Various alternatives were tried to predict end point carbon. Models 1-10, on the same 
data set (see Table 1 for summary of results ). The adequacy of models has been compared on 
the basis of standard deviation in Fig (17). It can be seen that in Models 1-10, the standard 
deviation has gradually decreased to approximately 8 milli percent. The independent 
variables finally selected for regression are time, lance height, ore (added as coolant after 
sublance measurement) and basicity of slag. From theoretical view point, all these variables 
do affect decrburization process. For example, lance height influences recirculation rate of 
metal, basicity of slag lying on top of metal also influences recirculation rate and residence 
time of droplets in slag and ore added directly supplies additional oxygen for decarburization. 
The respective signs of coefficients in Model No.’ 10 are also correct, positive for time, 
negative for lance height, positive for ore and negative for basicity. If average values of each 
of these variables is considered (time = 200 sec, lance height =2 m, ore added -300 kg, and 
basicity =3 ), then it would be seen for Model No, 10 that contribution of time - 1.72, 
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contribution of lance height=-0.015, contribution of ore =0.038 and contribution of basicity 
=-0.156. It shows that most important variables in converter operation is time, followed by 
basicity, ore and lance height. 

Comparative performance of carbon models of different plants (in terms of 
standard deviation of carbon) is shown in Fig (18) and the best performance is obtained for 
plant 2. 

Comparison of the performance of different temperature models (Modei No. 
15-24 summarized in Table {!)) is shown in Fig (19). The best results are obtained for 
Model No. 15 (ct = ± 4.986 "C) and Model No. 17 (cr = ± 4.986 '^C) for data set 1. The 
performance of same model when same variables applied to the data of different plants is 
compared in Fig (20); best results are obtained for plant 1. 

The variables selected in Model No. 15 and Model No. 17 are time (oxygen 
blown), lance life, temperature at sublance point, second blow additions, ore added (coolant), 
ratios of hot metal and scrap, end point carbon and basicity. From physical considerations 
also, all these variables are expected to affect bath temperature. Individual differences in the 
performance of models for different plants reflects upon the particular situation existing in a 
plant. 

To summarize : 

• the best model for carbon prediction for plant 1 is Model No. (10), a = ± 7.4204 milli %. 

• the best model for carbon prediction for plant 2 is Model No. (11), a = ± 6.8061 
milli %. 

• the best model for carbon prediction for plant 3 is Model No.(12), g = ± 7.6765 
milli %. 

• the best model for carbon prediction for plant 4 is Model No. (13), G — i 8.0459 milli % 
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• the best model for temperature prediction for plant 1 is Model No.(17), cr = ±4.986 °C 

• the best model for temperature prediction for plant 2 is Model No.(24), a = ± 5.518 °C 

• the best model for temperature prediction for plant 3 is Model No.(25), a = ± 6. 1215 °C 

• the best model for temperature prediction for plant 4 is Model No.(27), ct = ± 8.2032 

It may be noted that Model Nos. (9), (II), (12) and (14) for carbon prediction 
have been developed in the present work. These results show that carbon oxidation 
essentially occurs in the recirculation zone just below the jet. The best model for temperature 
prediction is a modified form of the model proposed by Bessho et al. and Byun et al. 



standard Deviation 
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Fig. 17. Comparative performance of diffrent carbon models (Model Nos. 1-10) of plant (1) 
Data set (compared values of y=mx curve, before ± 2o ) 







standard deviation 



Fig. 19. Comparative performance of different temperature models(Modei Nos. 15-23) of 
paint (1) Data set (compared values of y=mx+c curve, after ±2a ) 
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CHAPTER 4 

CONCLUSIONS AND SUGGESTIONS FOR FURTHER WORK 

4.0 Conclusions 

(1) . The models developed in I lie present work for prediction of end point carbon and 
temperature give better prcdiciioiis as compared to the exponential model and other models 
reported in the lileratiire. It is round that decarburization occurs primarily in recirculation 
zone in and around the area of iet impingement. 

(2) . Lowest possible standard deviations and heighest correlation coefficients, obtained at 

four different plants for the end point carbon prediction are 

Variables selected are 

Plant (1) a =±7.4204 time (0.0086379) 

R = 0.7728 lance height (-0.0076609) 

(for Model No. 10) ore added (0.0001 127) 

basicity (-0.052601 1) 

Constant (2.268365) 

Equation for prediction ofCr = C,. + (Csi, - Cc).exp{-[ Constant + time (0.0086379)+ lance 

height (-0.0076609)+ ore added (0.0001127) + basicity (- 


0 . 052601 1 ) 1 ) 
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Plant (2) a -± 6.S()6 1 time (0.008465324) 

R = 0.6509 ore added (0.000019422) 

(for Model No. 1 1 ) basicity (-0.0240933 16) 

Constant (0.490082) 

Equation for prediclion of Cr - i- (Csi, - ).exp{-[ Constant + time (0.008465324) +ore 

added (0.000019422) + basicity (-0.024093316)]} 


Plant (3) a -+7.6705 time (0.0 12579485) 

R =•- 0.66.1 1 ore added (0.000104090) 

(for Model No. 12) basicity (-0.035449106) 

Constant (0.403394) 

Equation for prediclion of C| = C'c + (Csi. -Cc).exp{-[ Constant + time (0.012579485) +ore 

added (0.000104090) + basicity (-0.035449106)]} 


Plant (4) a =±8.0459 time (0.008016913) 

R = 0.6472 ore added (0.00003738) 

(for Model No. 13) hot ratio (0.008102343) 

dolomite (-0.0001 12685) 

Constant (1. 122772) 

Equation for prediction of Cr =Cc +(C.si, -Cc).exp{ -[Constant + time (0.008016913) +ore 

added (0.00003738) + hot ratio (0.008102343) + dolomite (- 
O.OOOl 12685)]) 
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(3). Lowest possible staiiclard dcviaiion and heigliest correlation coettlcients, obtained for 
temperature prediction in the four plants are 


Plant (1) a “±d.9S() 

R - ().8d5S 


( Ibi' Motlcl No. I 5) 


Variables selected are 
lance life (0.0230456) 
temperature at sublance 
(0,6939638) 
oie added (-0.0093 125) 
basicity (0.291021) 
second blow additions 


(-0.0030052) 

hot ratio (1,2891021) 

A 02(7.436825) 

(1/Cr-l/CsL) (0.7156683) 
Constant (493.773952) 

Equation for prediction ot'Tr =('onslaiU + lance life (0.0230456) + temperature at sublance 

(0.6039638) +ore added (-0,0093125) +basicity (0.291021) 
■I second blow additions (-0,0030052) +hot ratio (1.2891021) + 


A (h (7.436825) + (1/Ci-1/Csl) (0.7156683) 


Plant (2) a --=±5.5 1 8 


R- 0.8785 


(for Model No. 24) 


temperature at sublance 

(0.86734136) 

ore added (-0,006905308) 


second blow additions 
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(-0.008240974) 
oxygen blown 
(0.036247155) 

Constant (210.817727) 

Equation for precliclion of 'ft -"Coiisianl -i- temperature at sublance (0.86734136) +ore added 

(-0.()()()905308) +second blow additions (-0.008240974) + 
o.xygen blown (0.036247155) 


Plant (3) a =-±6.12 1 5 temperature at sublance 

l< " 0.8712 (0.849737782) 

(for Mode! No.2.5) ore added (-0.004436813) 

second blow additions 


(-0.007871726) 
oxygen blown 
(0.037140187) 

Constant (234.027158) 

Equation for prediction ofTf --=Consiimt -t- temperature at sublance (0.849737782) +ore added 

(-0.0044368 1 3) -i-second blow additions (-0.007871726) + 
oxygen blown (0.037140187) 


Plant (4) a =± 8.2032 

R = 0.8575 


temperature at sublance 
(0.881824165) 


(for Model No.27) 


ore added (-0.010025682) 



100 


lime(-0.000257262) 
dolomite (-0.02135692) 

A 02(12,078419703) 

(CsL-Cr) (-0.040002643) 

Constant =221.672589 

Equation for prediction of'lV “Constant + temperature at sublance (0.881824165) +ore added 

(-0.0 10025682) -i-lime (-0,000257262)+dolomite (-0.02135692) 
t A (>2 (12.07X4 1 970.1) -I- (C| -Csi.)(-0.040()02643) 

(4) . Both lance lieiglil and lance lil'c arc important variables and if data on these are made 
available then carbon prediction can be improved by nearly 0.8 milli percent and temperature 
prediction by 0.8 "C. 

(5) . Slag basicity is a more significant variable than slag volume. The former determines 

the viscosity and hence foaming characteristics of slag. 

(6) . Slightly better prediction of temperature and carbon can be obtained in converter 

equipped with bottom stirring (plants 1.2 and 3 ) than otherwise (plant 4). 

4.1 Suggestions for Furtlier work 

(1). 'I'he models developed in the present wtirk need to be extended to the prediction of 
end point manganese and phosphorus as well. Preliminary work has already shown that 
phosphorus can be predicted in the range of o 1.6 milli percent and manganese in the 
range of a=± 5 milli perccnl. 



(2). Specific inlluciiccs oflnncc height anti noitzle design on process kinetics need to be 
studied in detail. 


(3). The inodchs obtainetl in the present work can be Hii'ther tuned by using soft computing 
techniques including neural net:;, genetic algorithms and fuzzy logic. 


Table. (I) Siiniinnry Dfic.Miihs nTnuulcIs for caiboii aiul tcinpcrnluie prediction 
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Mo 

del 

No. 

Dependent 
variable in 
regression 

Indcpendcnl 
variables taken 
in regression 

Selected variables by 
stepwise regression and 
corresponding cocflicients 

For the 

regression 

curve 

Before ±2a 

After ±2 a 




Rc.su lt.s for E.xponential 
model (plant (1 ) data set) 




1 

Cp 

teinpenttinc at 
sublancc, 
second blow 
additions, 
slag volume, 
ore added, 
hotratio. 

orciicldod(-().()tlO()5.12), 

holnilio((),01«72(m), 

Regression constaiil 
(0.-tH4(my). 

y=m,\ 

11=301 
111=0.98933 
a =14.002 
R=0.66119 
F=6865 

11=284 
111=1.00528 
a = 11.4872 
R=().73016 
F=9404 

y=mx+c 

n=30l 

111=0.5981 

c=27.3968 

ft = 12.0415 

R=().6612 

F=232 

11=285 

111=0.5958 
c=27.002 
a =10.1623 
R=0.7l 

F=287 

2 

Cp 

lance height, 
tcmperatuie at 
sublance, 
second blow 
additions, 
slag volume, 
ore added, 
hotratio 

liincc Iieiglit(-I).IH)()27K5), 

orcaddcd(-0.(l0(H)372). 

lmliiitio((),()lK5K(i()). 

Kgrcssioii constaiil 
(0.509299). 

y=m.\ 

11=301 
111=0.98808 
a =13.7495 
R=0.66998 
F=7t02 

11=292 

111=0.99813 

ct«12.1226 

R=().7l()26 

F=8724 

y=m.x+c 

11=30 1 
111=0.5782 
c =27.3066 
a =11.7602 
R=0.67 
F=243 

11=287 

111=0.5809 
c =28.0326 
rr =9.9284 
R=0.708l 

F=286 

3 

Cp 

All the variables 
(temperature at 
sublancc. 
second blow 
additions, 
slag volume, 
ore added, 
hotratio) less 
thier 

corresponding 
mean values. 

(ore added-mcan of ore 
added) 

(-().()()(){)5.12), 

(hotrntio-mcan of 
hotratio)({).OIH72()H), 
Rcgrc.sslon cons(ant(0.55,S6). 

y=m.v 

11=301 

111=0.96853 

fT=13.8538 

R=().66119 

F=6865 

n=286 
111=0.95147 
a = 11.2682 
R=0.7349 
F=9373 

y=mx+c 

11=301 

111=0.7309 

c= 16.6945 

or=13.31 

R=0.6612 

F«232 

11=284 

111=0.69 16 
c= 17.6208 
a =10.3263 
R=().7301 

F=321 




Sublancc model iMethod 1 
(for plant ( 1 ) data set ) 




1 

K 

temperature at 
sublancc, 
second blow 
additions.ore 
added, hotratio 

ore ndded(0,b'oooi)iif).f5), 
biisieil,Y(-0.00()2 15759). 
Kegrc.s.sion consinni 
(0.10422). 

y=mx 

n=301 
111=0.9542 
a =13.33 
11=0.6627 
F=7046 

n=28H 

111=0.96759 

a=n.43 

R=().73l5 

F=9130 

5 

K 

lance height, 
temperature at 
snblancc,.second 
blow additions, 
ore 

added.hotratio 

Iniice liciglil(0.li()00().54F)7), 
ore:iddod(0.000000(.K.)). 
bnsiei(y(-0.()002l9421). 
Uegresion con.stnni 
(0.0101.55). 

y=mx 

11=301 
111=0.9539 
a =13.173 
R=0.6688 
F=72 1 1 



* n=nuinbcr of (lain points. ni“.slo|)c orciin c. c=constanl. ^.standarc! dcviatioiu !l-corrclalion 
!cocnicicnl.P=F-ic.s( value 


(continued ...) 
































Summary of results of models for carbon and temperature prediction (continued ..) 


10 

C -C 

In L * - 

timcdance 

heighl, 

temperature at 
siiblancc.second 
blow additions, 
ore added, 
hot ratio.bnsicity 

time{0.()()86379), 
lancc height(-(). 0076609), 
ore added(0.000 1127), 
bnsicity(-0.05260||). 
Regression constant 
(2.26H365). 

y=mx 

11=272 
111=0.94839 
a =11,61 
R=0.74384 
F=8341 ' 

n=263 
m=0.96813 
a =10.0727 
R=0.78413 

F= 10739 

111 — 

y=mx+c 

11=272 

111=0.5170 

C“30.3324 

a =8,3772 

R=0.7438 

F«334 

n=262 

111=0,5284 

c-29.3979 

CT =7.4204 

R=0.7728 

F«385 




Modified method (2) 

(for plant (2) data set ) 




11 


time,tcmperaturc 
at sublancc, 
second blow 
additions, 
ore added, 
hot ratio.bnsicity 

liiiie(0,(mH465:?24). 
ore !idclcd(0.0(l()0 19422). 
b;isicily(-0,()24t)933 16). 
Regression coiisliiiil 
(0.49()OK2). 

y=mx 

11=1744 

111=0.96204 

0=10.2214 

R=0.6341 

F=53258 

11=1667 

111=0.98114 

CT =8.7025 

R«0.6642 

F«698I8 

c - c 

y=mx+c 

11=1744 
111=0.4664 
c=29.9053 . 
rr =7.6145 
R=0.6341 
F=1171 

n=167l 

111=0.4545 

c=30.5902 

CT =6.8061 

R=0.6509 

F=1227 




Modified method (2) 

(for plant (3) data set ) 




12 

C -C 

In r - - 

time.tcmpcralurc 
at sublancc, 
second blow 
additions, 
ore added, 

Itot raiio.basicily 

lime(().OI2.‘>794K.S). 
oreadded(().OOOI()4l)9()). 
basicily(-0.(l35449106). 
Regression coasinni 
(0.403394). 

y=mx 

11=748 

111=0.9491 

CT =12.2767 1 

R=0.6192l 

F= 16771 

11=715 

111=0.97358 

CT =10.0979 
R=0.6363l 
F=23409 

in 

C -C 

y=mx+c 

11=748 

111=0.4343 

c=32.6684 

a =8.8244 

R=0.6192 

F=463 

11=711 

111=0.4329 

c=32.4388 

CT -7.6765 
R-0.663 1 
F=556 


* n»mnnbcr ordain polnis. in=slo|)C of cnn'c, c~conslnnl, “lanndnrd dcvlalion, R=cofrclnlioii (continued ...) 
co£rrictcnlJ"=F-icsl value 
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Summary of results of models for carbon and temperature prediction (continued ..) 


6 

K j 

1 

) 

i 

i 

1 

lance height, 
temperature at 
sublaiicc.sccond 
blow additions, 
slag voluinc.orc 
addcd.holratio 

lance height 

(0.0000032X8). 

ore addcd(0.000000764). 

Regression constant 

(0.009085). 

y=m.\ 

n=10l 
m=0.95329 
a -13.522 
R=t).6497 
F-6834 

! 


1 

Snblancc model: Method 2 
(for plant (1) data set ) 



1 

i 

7 

'• tiinc.icnipcratiirc 
j at siiblance. 

C * second 

— In blow additions. 

^ SL ‘ orcadded. 

! hotralio. 

1 basicity 

111110(0.007827967). 
ore addcch 0.000 1 14356). 
basicilv(-0.04003613 1 ). 
Rcurc.ssion constant 
(0 473776). 



y=in,\ 

11=30 1 
tn=0.953ll 
rT= 12.9347 
R=0. 66582 
F-7466 

n=289 1 

ni=0.9T334 j 
a =11.027 ! 
R-0.73165 i 
F=9867 i 

>-m.\+c 

n=.>01 

m=t).496 

c-32.0314 

a = 8615 

R-0.6658 

F-238 

n=282 

m=0.497 : 

c=3l.4l41 j 
rr =8.3363 | 

R-0.7138 1 

F=290 ! 

8 

tinic.lancc 

height. 

C tcinpcniiurc at 

In sublaiicc.sccond 

^'SL blow additions. 

: ore added, 

, hot ratio.basicity 

I 

(itnc(0.0078753), 
lance hciglit(4).0()66976), 
ore addcd(().000 1102). 
basicity(-0. 0466949). 
Regression constant 
(2.013246). 

i 

y=m.\ 

n=3()l 
in=0, 95856 
rr «11.816 
R-0.73992 
F=8226 

n=26 1 i 

in=0. 98009 1 

0=10.012 i 

R=0,79227 1 
F= 11022 j 

N-m.\"K: 

n=272 

m=0.5692 

c=27.3729 

n = 9.333 

R=0.7399 

F=326 

n=259 i 

m=().5977 

0=25.5638 

a = 8.092 

R-0.7784 

F=395 


j 

Modified method 2 
(for plant (1) data set ) 




9 

^ SL ^ If 

time, 

temperauirc at 
sublancc.sccond 
blow additions, 
ore added, 
hot ratio.basicity 

timc((l.(l()K5.')4l7()). 
oic:iddccl(tl..0()<»IIR.'i46). 
bnsicjly(-0.()4342846). 
Regression constant 
(0.505659). 

y»m.\ 

11=30 1 
m=l).94268 
cr =12.6652 
R=0.66826 
F-7618 

n=286 
m=0,967()2 
a =10.5013 
R=0.73212 
F=10511 

\'=mx+c 

n=301 

m«0.4428 

c-35.0092 

(T =8.7458 

R-0.6683 

F=241 

n=286 

m=0,45 

c=34.1851 

a =7.613 

R=0.712 

F=29i 


— ^ ^continued'...) 

* n=iuuiibcr of data points. ni=slopc ofcim e. c=coi!StaiU, =siandard deviation, R-correlation 
cocfricicnl.F=F-tcst value 
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Summary oF results of models Ibr carbon and temperature prediction (continued .,) 





Modillod method (2) 

(for plant (4) data set ) 




13 

-In 

timc.tempcrnturc 
at sublancc, 
hot ratio.orc 
added. basicity, 
dolomite 
(end point 
carbon<().06) 

timc(().0()St)U)'>i3). 
ore iHlc!e«l(0.00(ll)374K), 
luilmlio((l.(l()K 102.44.1). 
dc)loniiic(-(>.0001 laoK.I), 
le^'res'sioii coii.siiiiil 
(1.122772). 

y=m,\ 

11=1832 
111=1.01352 
a =9.5153 
R=0.58892 
F=122I6 

n=1742 
111=1.01989 
a =8. 1368 
R=().6739 
F=443I6 

y-m.N+c 

11=1832 

m=0.75l8 

0=10.8853 

a =9.21 16 

R=().5889 

F=97I 

11=1755 
111=0.772 
c= 10.0659 
a =8.0459 
R=0.6472 
F=1263 

14 

C -(’ 
-In ' " 

time, temperature 
at sublancc, 
ore added, 
hot ratio. basicity, 
(lime+doloniile) 

timc(0.0()S209.S46). 
basicilyto. 024547185), 
(lintet dolomite) 
(-O.()0OOU502). 
hotralio(0.(K)7K97245), 
Regression constant 
(0.990858). 

y=mx 

11=18.32 
m= 1.0201 
a =9.475 
R=(),6I437 
F=34958 

11=1750 
m= 1.02497 

CT =8.1837 
R=0.69078 
F=44558 

y=mx*fc 

11=1832 
111=0.8092 
0=8.7711) 
a =9.279 
R=0.6!44 
F=no9 

11=1758 

m=().8233 

c=H.ll64 

a =8. 1364 

R=0.667 

F=1407 




Tcniponilurc model 
(for plant (l)dala set ) 




15 

End point 
temperature 

lance life, 
tcmpcnilurc at 
sublance..sccond 
blow additions, 
ore added, 
hot ratio, basicity, 
AO. 

(l/Cr-l/Csi.) 

lance lire(().02.10456), 

temperainre at 

.siihlaiicc(0.(>9396.1H). 

.second blow additions 

(-0.00.30052), 

ore addecK -0.0093 125), 

hot ralio( 1.289 1021), 

AO:.(7.436825), 

(I/Ct-1/C.si.)(0.7l56683), 

basicity(0.29l02l), 

regression con.siant 

(49.^.773952). 

y=mx 

11=301 
111=0.99998 
a -7.4481 
R=0.79H88 

F= 15068846 

11=284 
111=0.99969 
a -6. 1353 
R=0.857 
F=2()943815 

y=mx+c 

11=301 
111=0.63816 
c=602. 99247 
nr =5.9596 
R=0.798H8 
F=527 

n=287 

111=0.6442 

c=592.4415 

a =4.986 

R=0,8458 

F=716 

16 

End point 
temperature 

lance height, 
tempcralurc at 
sublance,.second 
blow adtlition.s, 
ore added, 
hot ratio.slag 
volume, AOj 

(l/Cf-l/C'st. ) 

temperature at 
.sublancc(0. 708281 1), 
.second blow additions 
(-0.00.39HI 1), 
ore addod(-0.0()877.37), 
AO. (7.7 104538), 
(l/Cr-l/(b(.)(0.746I434). 
rcgrcs.sion constant 
(477.550594). 

y=mx 

11=301 
111=0.9998 
a =7.6108 
R=().794()7 

F= 130.37529 

11=256 
111=0,9997 
a “6. 1632 
R=(),85716 

F= 18704929 

y=inx-t-c 

11=301 
m“().63054 
c=6 15.6 147 
a =6.0546 
R=().79407 
F=46t) 

n=258 
111=0.6250 
c *624 
a =4.8919 
R=0.8458 
F=643 


’•ii»nninbcr of data points. in«slopc ofairvo, c-constanl. -standard deviation, R- nrclation 
ci;cfricicnl,F*F-tcsl value 


(continued ...) 




106 


Summary ofresulls of models tor carbon and temperature prediction (continued .,) 


17 

End point 
Icinpcrature 

lance 

hcighiJancc 

lirc,iempcratiire 

at 

sublancc.sccond 
blow additions, 
ore added, 
hot ratio, basicity, 
AO 2 

(1/Cr-l/Csi.) 

lance lile((r()2.UM5(i), 

(emporature at 
sublancc(0.()y3%3S), 
second blow additions 
(-().()030()52), 
oreaddcd(-n.()Oy3l25), 
hot raiio(1.2K‘> 11)21), 
A02(7.43bK25), 

(I/C, -1/C.si.)(t).715r)683), 
basicity(0.291()2l). 
regression constant 
(493.773952). 

y~- in.K 

n«30l 

111=0.99998’ 

CT =7.448 1 
R=0.79888 

F= 15068846 

n-2H4 
m=0. 99969 
ct=6.I355 
R=0.857 
F=2()943815 

y=m.\+c 

n=30l 

in=0.638l6 

c^6i)2.99241 

rr =5.9596 

R=0.79888 

F=527 

n=287 

m=0.6442 

c=592.44l5 

a =4.986 

R=0.8458 

F=716 

18 

End point 

lancc 

hu\ce lire(t).l)24ri530). 

y=mx 

n=3()l 

11=285 


lempcnUuie 

height, lance 

icinpcrature at 


m=0.99998 

m=0.99966 



lifc,iempcraturc 

sublancc(().7n5224). 


a =7.9 11 

a =6.68 



at 

second blow additions 


R=0.76931 

R=0.82328 



sublancc.sccond 



F= 13356695 

F= 17727552 



blow additions. 

oiciKklal(-(M)(mf.5f.l). 






ore added. 

A():,(y,2 189470), 

y~m\4c 

n=301 

n=284 



hot ratio, hasicily, 

icgrcsion constant 


111=0.5918 , 

m*0.59<)0 



AOj 

(47I.K75.‘i(m). 


c=r)K0.2434 

c=667.7K23 



(l/Csi.) 



a =6.0959 

CT =5.0482 






R*0.7693 

R=0.8272 






F=433 

F=611 

19 

End point 

temperature at 

tctiipcratiirc at 

y=m.K 

n=301 

u=285 


temperature 

sublancc.sccond 

siiblaiicc(().7 13 1490), 


111=0.99998 

111=0.99965 



blow additions, 

second blow additions 


a =7.997 

a =6.7517 



ore added. 

(-0.0037001 ). 


R=0.7635 

R=().82009 



hot ratio.slng 

ore added(-0.009()445). 


F= 1307 1322 

F= 1760 1267 



volume, AO 2 

AO:.(9.30()2509), 






(1/Csi.) 

rcgrs.sion constant 

y=mx+c 

n=301 

n=287 




(470.%7HKI). 


111=0.5829 ' 

m=0.5883 






c=695.0989 

c=685.612() 






CT =6.1 

(T=5.152 






R=0.7635 

R=0,8I67 






F«417 

F=570 

20 

End point 

lancciilc.lcmpcra 

1 a n c c i i i c ( ( i , i ) 2 • i f > 5 1 0 ) , 

y^’in.s 

ir-.ldi 

11=285 


temperature 

lure at sublancc, 

Icmpeialuic at 


111=0.99998 

111=0.99966 



second blow 

.sublaMce(0.7 115224). 


(T =7,911 

CT =6.68 



additions. 

second blow addilioiKS 


14=0.7693 1 

R=0.82328 



ore added. 

(-O.ooni-I.S), 


f'= 13.3 56695 

F= 17727552 



hot ratio.slng 

ore addcd(-0. 008050 1). 






voluiuc, AO-. 

A():.{9.2IH9470). 

y=inx+c 

n=30l 

11=284 



(l/Cs,.) 

regresion constant 


111=0.5918 

111=0.5990 




(47I.H7550H), 


0=680.2434 

0=667.7823 






rr =6,0959 

a =5.0482 






R=0.7693 

R=0.8272 



1 



F=433 

F=6ll 


■* n*n»nnbcr of dnln points. in*slopc of curve, c”constant, ^standard deviation, ll=canclalioii 
•cocfricicnl.F®F-!csl value 


(continued ...) 








Suminary of results of models for carbon and temperature prediction (continued ,.) 


21 

End point 
temperature 




lance lieiglit. 
temperature at 
sublaiice.sccond 
blow additions, 
ore added, 
hot nilio.slag 
volume, AO: 

(1/CsO 

lance lieiglit(-0. l(M959;i), 
temperature at 
sublance(0. 7 1390X1), 
second blow a(lditj{)ns 
(-0.t)()456(»5), 
ore addcd(-().OOX3f)9). 
A():(‘).395H17()). 
regrcs.sioii con.staiU 
(492.736455). 

jj 

n=30l 
m=().9998 
a =7.9949 
R=0.76962 

F« 11814652 

11=257 
111=0,99964 
a =6.6299 
R=0.83I8() 

F“ 16226661 

y=inx'ic 

n=3()l 
10=0.59231 
c =679.3 14 
a =6. 1645 
R=0.76962 
F=392 

11=258 

111=0.6084 
c =651.8922 
0=5.1054 
R=0.8301 

F=567 

22 

End point 
temperature 

icmpcralurc at 
sublnncc.second 
blow additions, 
ore added, 
hot raiio.slag 
volume, oxygen 
flow rale, 

(l/Cm.) 

icmpcraturc ni 

sublance(0.()9Xl714). 

second blow additions 

(-O.Oon.SoX). 

ore added(-0.0092944). 

hot ratiot 1.6559 1()7), 

oxygen How rale 

(0.0276419). 

regre.s.sion 

conslanl(4HH. 152365) 

y=mx 

n=30l 

111=0.99998 

CT =8.0494 
R=0.759H9 

F= 1290 1493 


23 

End point 
temperature 

lance life. 
icmpcnUurc at 
sublance.sccond 
blow additions, 
ore added, 

Iiot ratio. slag 
volume. AO^ 

(1/Cf-l/Cs,.) 

lend point 
carbon 

calculated from 
carbon model 
No.(K) 

lancc Iife(0.()246530), 
temperature at 
sublance(0.7l 15224), 
second blow additions 
(-0.0037345). 
ore added(-0.00H656l), 
A():.(9.2 1X9470). 
rcgrcsioii constant 
(47I.X755()X). 

y=mx 

11=301 
111=0.99998 
a =7.911 
R=0.7693I 

F= 13 356695 

n=285 

111=0.99966 

0=6.68 

R=().82328 

F= 17727552 

y=m.\4-c 

11=301 

111=0.5918 

c=680.2434 

a =6.0959 

R=(),7693 

F=433 

11=284 

111=0.5990 

c=667.7823 

CT =5.0482 
R=0.8272 

F=6I1 




Tcmparaluc model 
(for plant (2) data det ) 




24 

End point 
temperature 

tcmpcniturc at 
sublance.sccond 
blow additions, 
ore added, 
hot miio.slag 
volume, oxygen 
How rale. 

(I/ChiJ 

(cmpcratiirc at 
sublance(0.X673413r)), 
second blow additions 
(-0.00X240974). 
ore nddcd(-0.()()69()530X), 
oxygon How rate 
(0.03(1247155). 
regre.s.sion con.sljml 
(2 10. XI 7727). 

y=mx 

11=1744 
111=0.99998 
a =7.5264 
R=().84H34 
F=86()61195 

11=1659 

in-0.99996 

CT =6.5335 
R=().H7849 

F= 108632432 

y=m,X'l-c 

n=1744 
m=0.7l97 ' 
c=468.72()3 
n = 6.3866 
R=O.H4X3 
F=4472 

n=1662 

m«0.7I29 

c=48{),0948 

CT =5.518 
R=0.H785 
F=5f)H 


’ n»mimbcr of dnin poiiKs. in*slo|)c of ctin'c, c~consl;uii, ^slaiulnrtl dcvialion, R==corrclalion 
coefficient. F~F-(cs( value 


(continued ...) 
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Summary ol rcsuKs ornuKlels for carbon and temperature prediction (continued ..) 



Htici point 
Icinpciotnrc 

Icinpcintinc at 
snblancc. second 
t)lo\v luldilions. 
ore ndded. 
hot intio.slofA 
\ (dnnic. oxN^cn 
IUn\ lolc. 

(l/('si) 

'I'empetalmo niotlcl 

(Ibr plant ( 0 data set ) 

letnperalnie al 

.stil)lance(().KI')7n7K2). 

second tilow additions 

(-l),oi)7X7177.r>). 

ore added(>0 ()l)447f)8Ll). 

oxvp.en flow rale 

(0.0,17141)1X7). 

lemession constant 

(2.I4,027LSX). 




25 

>'='in.\ 

n=748 

m=l).99997 

rr =8.613.1 

R=0.829()8 

F-28262053 

11=710 

111=0.99999 
a -7.2951 
R-0.87052 
F=37386359 

y=mN+c 

n=748 

m=().6R73 

c=523.5034 

(7=7.1454 

R=0.8291 

F=1640 

11=712 

111=0.7041 

c=495.H79 

<7*6.1215 

R=0.8712 

F=2235 




'remperatme model 







(for plant (4) daht set ) 
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End point 

tcinpcratnrc at 

ternperatnre at 

>’=mx 

n=22i4 

11=2100 


tempentture 

siiblancc, 

Mil)laiia'((),8(»2(>()4,t Ki). 


m=0. 99992 

111=0.99968 



ore added. 

ore added! -0.0 KWyTW)). 


a =12.0.45 

0=10.1074 



(I/C'f - 1/t'si,). 

(l/t'i - l/Csi.) 


R=0.8()6()2 

R=0.85358 



lime. 

(1-^07. 1725 1.5070). 


F=41035957 

F=55202()6l 



dolomite. 

Iinic(-0,(){)4517I0). 






(Csi.-C'i) . AO: 

((‘.!,.-(‘,)(-0.02292:I%H), 

y^mx+c 

11=2214 

11=2111 




AO: (H.5b 1022479). 


111=0.6498 

111=0.6607 




regression conslanl 


c=.574.7695 

C=555,6642 




(251.K2195I2). 


a =9.70.52 

a *8.2886 






R=l),X060 

R=0.8477 






F=4lOI ' 

F=5384 

27 

End point 

temperature at 

tcinpcralmc al 

y=mx 

11=1819 

11-1723 


Icinpcratnrc 

sub la nee. 

snblancc(O.XKIH24t65), 


111=0.99993 

111=0.99966 



ore added. 

ore addcd(-0.0 10025682), 


nr = 11.9799 

a =9.8743 



(l/Cr - 1/Csu). 

(1/C| - l/fsi.) 


R=0.8033 

R={).85616 



lime. 

(-7fi(..2.‘'lX5IKOl), 


F=3 4003698 

F=47425417 



dolomite. 

limc(-0.0002.';72(.2). 






(('si.-('i). AO: 

(('ia -('d (-0040002(.4.7). 

y=mx+c 

11=1819 

11=1738 



(cat bon 

A().>( 12,07X4 10704). 


111=0.6452 

m=0.6685 



calculated Irom 

doluinile(-0.02 1.45002). 


c=58 1.0997 

c=542.7269 



carbon model 

rcuression coiisiniu 


a =9.6253 

a =8.2032 



No.(14) 

(221.6725X0), 


R=0,8033 

R=().85l5 






F=3305 

F=4579 


' imtmmbcr nCdnla poiMl<;, iir slope orciiive. c^coii.slanl. -sliiiulaid dcvialion. lt=coi relation 
cocllieicnl.l-' r-K'sl value 
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APPENDrX-1 


Various operational variables given in tlie sample data sets are calculated as follows. 


Time of blow can be calculated as described inChapter3 byEqs. (3.11) or (3.12). 


Baxicitv 


Cap 

Sip. 


( 'up 


StUn.HM )xWl.()f.HMx — . - 

28 100 


AO, = 


0„ 


ij'Vt. of .HMxO. 92 + IVf. of .Scrap) 


HutRatio = 


Wl.of.HM 
Wt. of .Scrap 


Slagvol = limel+Rslagl+dolol+dolspI I+Si02 

Second blow additions = lime2 + Rslag2 +dolo2 -r-dolo2+dolspl2 

Abbreviations used are : 

HM =hotmetal 

SBA = second blow additions 
limel =lime added in the first blow 


lirne2 = lime added in the second blow 



dolol =dolomite added in the first blow 

dolo2 = dolomite added in the second blow 

dolspll = dolomite special added in the first blow 

dolspl2 = dolomite special added in the second blow 

Rslagl = return slag in yhe first blow 

Rslag2 = return slag in the secon blow 

Orel =ore added in the first blow 

ore2 = ore added in the second blow 

O;: = oxygen blown in the second blow 



sitiiiplt' (liU.ii sot of i)laiili(l) 


Hoat, NoT] 

liolniotal 

si(hriiM)T 

»S(na|) 1 
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t7 T 
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tons 

'mini %) 
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®C < 
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4' 

9nilli “/ 
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1636 ^ 
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61 
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70 
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57.7 

2420 

1593 
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67 
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228 
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1623 
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60 
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64 
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218G 
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75 
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482 
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2H)2 

279.0 
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2000 
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502 
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62 
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300.0 
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1584 
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67 
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61.4 

1950 
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451 
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62 

8510 
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296.6 
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54.6 

2340 

1612 

472 
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54 
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2199 
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1581 
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49 
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13070 

2213 

285.0 
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1700 
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1672 

59 
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460 

55.5 
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1602 

356 

1648 
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441 

60.4 

1680 

1613 

358 

1672 

69 

10800 

2218 

285.4 

447 

53.7 

1980 

1621 

386 

1682 

56 

10230 

2219 

285.3 

408 

46.0 

1620 

1607 

442 

1663 

79 

9620 

2220 

285.0 

384 

57.9 

1440 

1614 

285 

1650 

97 

11250 

2221 

285.0 

348 

58.0 

1770 

1614 

377 

1679 

71 

10120 

2223 

285.2 

636 

44.5 

1680 

1599 

234 

1668 

45 

14220 

2224 

285.0 

' 389 

59.3 

1280 

1620 

251 

1653 

87 

11630 

2225 

285.0 

447 

63.5 

1300 

1622 

354 

1661 

85 

11020 

2227 

275.3 

485 

57.7 

1500 

1610 

411 

1657 

82 

11100 

2228 

281.3 

454 

58.5 

1610 

1624 

409 

1664 

119 

13060 

2229 

284.1 

548 

66.6 

1320 

1628 

351 

1666 

76 

12490 

2230 

283.8 

671 

72.3 

1100 

1608 

292 

1643 

100 

14800 

2231 

282.2 

383 

58.5 

1980 

1595 

470 

1638 

77 

9300 

2232 

280.2 

375 

63.8 

1990 

1580 

442 

1658 

65 

9210 

2233 

292.8 

358 

47.6 

2100 

1627 

550 

1671 

88 

10260 
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IJcat, No. ,Sla.gvol I lime 


1 5826.9082 
14343.3145 
13672.2861 

14534.1719 
13573.7422 

17478.8574 

18381.8574 
14968.8916 
1 7759.9707 

15337.1719 
14805.4287 
14317.4854 
16089.7012 
13360.0488 
15707.2227 
15734.6426 

18180.2852 

16153.8574 
19136.0352 
20189.9277 

22209.8926 
19400.0723 
19578.8457 
17741.4277 
19929.5430 
16604. .5898 
1.5967.7246 
16434.3359 
16409.1426 

16345.2852 
22226.8691 
17307.6777 

16969.8926 
17161.1.543 
18800.6465 
19446.1465 
23180.6387 
1.5986.0557 
16.551.6074 
16014.1943 


3.7733 

4.9120 

3.9264 

4.5195 

3.9932 

3.7421 

3.8567 

4.7971 

4.7617 

3.6585 

4.8954 

4.0368 

4,7723 

3.7438 

3.6269 

4.0154 

4.0904 

4.5677 

4.10.52 

4.2537 

4.4084 

5.1236 


6.0709 

3.7431 

4.4545 

5.8069 

3.9923 

5.2115 

4.8947 

5.1899 

6.2605 

5.2516 

4.4986 

5.5278 

5.4744 

3.9813 

3.9914 

4.8235 

5.0741 

4.0247 

3.2994 

6.2240 

6.1881 

6.6251 


13.7825 
5.1028 
14.2825 
6.8064 
11.6995 
15.2665 
10.3958 
6.8005 
11.4320 
17.9487 
7.5102 
8.9398 
5.9584 
10.3089 
6.5753 
10.8594 
13.1222 
9.. 54. 55 
9.8457 
14.6707 
8.7245 
10.1958 
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APPENI)lX-2 


Sample Data set of plant (2) 


(for details of variables see Appendix 1) 


Heat No. 

Ilotinetal 

(kg) 

Sii 1 
[in HM] 
,% 

Scrap 

(kg) 

Cs>. 

% 

0 22 
(nv') 

Ore2 

(kg) 

11127 

322000 

0.205 

9704 

0.210 

1692 

0 

11153 

323100 

0.306 

10921 

0.222 

1742 

0 

11165 

325300 

0.326 

9798 

0.347 

1951 

6 

111()7 

322700 

0.270 

9430 

0.388 

2520 

0 

11186 

319700 

0.205 

1 1.534 

0.520 

2339 

0 

11202 

330400 

0.320 

9497 

0.478 

2260 

0 

1 1213 

314100 

0.438 

14483 

0.362 

1618 

0 

\\2\4 

330900 

0.507 

14016 

0.339 

1545 

0 

11264 

331100 

0.342 

13695 

0.391 

1531 

0 

11282 

334200 

0.266 

9488. 

■0..328 

2038 

0 

1 1305 

316500 

0.356 

10147 

0.310 

1697 

0 

11306 

330500 

0.272 

10104 

0.431 

2027 

0 

1 1 3 1 1 

325800 

0.382 

10208 

0..506 

1995 

0 

11326 

325400 

0.240 

10985 

0.633 

2225 

0 

11330 

319800 

0.667 

13601 

0.325 

1972 

0 

11331 

333100 

0.246 

12027 

0.349 

1921 

0 

11332 

330900 

0.266 

10992 

0.265 

1824 

0 

11351 

332100 

0.315 

11998 

0.269 

1718 

0 

11366 

328200 

0.567 

13315 

0.428 

2062 

0 

11371 

336000 

0.283 

9122 

0.344 

1930 

0 

11386 

313700 

0.241 

9711 

0.393 

1947 

0 

11389 

316200 

0.647 

15242 

0.416 

1733 

0 

11398 

330900 

0.372 

128.55 

0..397 

2137 

0 

1 1425 

314500 

0.279 

10169 

0.354 

21,53 

0 

11426 

315800 

0.318 

13743 

0.494 

2180 

0 

1 1427 

323400 

0.472 

1466.5 

0.360 

2139 

0 

11490 

330700 

0.302 

12200 

0.290 

1833 

0 

11511 

329300 

0.336 

11931 

0.446 

2030 

0 

11533 

331800 

0.300 

9397 

0.309 

1676 

0 

11534 

335100 

0.350 

1 1243 

0.390 

2012 

0 

11535 

316600 

0.303 

9595 

0.289 

1522 

0 

11555 

325200 

0.305 

9.591 

0.355 

1821 

0 

11556 

321700 

0.442 

11684 

0.287 

1783 

0 

11565 

321000 

0.478 

13.539 

0.539 

2339 

dQ 

11569 

320400 

0.448 

11472 

0.405 

2136 

0 

11574 

328500 

0.329 

1 1646 

0.431 

2189 

0 

20019 

326600 

0.277 

9.591 

0.288 

1417 

0 

20056 

325200 

0.233 

10960 

0.436 

2257 

0 

20095 

323300 

0.316 

9495 

0..340 

1830 

0 

20096 

325200 

0.337 

9.504 

0.326 

1827 

0 


(coiUinucd ...) 





Sample Data set of plant 0.) (tontiiuied ...) 




time 

(sec) 


Tr 

("O 

SBA 

(kg) 

' Tsi. 

("C) 

11127 

0.061 

169.2 

6.8604 

1662 

0 

1602 

1115:5 

0.035 

174.2 

5.1548 

1661 

0 

1600 

11 1()5 

0.060 

195.1 

4.3116 

1683 

0 

1613 

111()7 

0.044 

252.0 

5.0507 

1676 

0 

1571 

11186 

0.057 

233.9 

8.2128 

1664 

0 

1575 

1 1 202 

0.052 

226.0 

4.1918 

1671 

7 

1595 

11213 

0.075 

161.8 

4.9127 

1654 

0 

1591 

112H 

0.099 

154.5 

3.8988 

1648 

0 

1599 

1126d 

0.086 

153.1 

5.6440 

1660 

0 

1611 

11282 

0.049 

203.8 

4.9807 

1683 

0 

1602 

11305 

0.060 

169.7 

4.2026 

1667 

0 

1611 

11306 

0.048 

202.7 

5.24.52 

1702 

0 

1614 

11311 

0.082 

199.5 

3.8277 

1651 

0 

1578 

11326 

0.100 

222.5 

6.5641 

1675 

0 

1594 

11330 

0.057 

197.2 

2.9756 

1692 

0 

1612 

1 l:53l 

0.052 

192.1 

6.8494 

1670 

62 

1.589 

11332 

0.063 

182.4 

5.8278 

1662 

2 

1600 

11351 

0.072 

171.8 

5.3522 

1663 

3 

1605 

11366 

0.054 

206.2 

3.3.391 

1670 

0 

1595 

11371 

0.050 

193.0 

4.4768 

1677 

0 

1599 

11:386 

0.080 

194.7 

5.9943 

1666 

0 

1599 

11389 

0.072 

173.3 

3.4768 

1668 

0 

1610 

11398 

0.047 

213.7 

4.8735 

1675 

0 

1599 

11425 

0.049 

215.3 

5.4083 

1656 

0 

1578 

11426 

0.082 

218.0 

6.3863 

1669 

0 

1598 

11427 

0.058 

213.9 

4.4834 

1671 

0 

1597 

11490 

0.056 

183.3 

5.7007 

1678 

4 

1622 

11511 

0.072 

203.0 

5.0321 

1653 

0 

1570 

1 15:53 

0.052 

167.6 

4.40.55 

1675 

0 

1602 

11534 

0.053 

201.2 

4.4735 

1699 

0 

1609 

1 1535 

0.058 

152.2 

4.6676 

1663 

0 

1615 

1 1555 

0.048 

182.1 

4.5125 

1659 

0 

1592 

1 1556 

0.042 

178.3 

3.8346 

1670 

0 

1614 

11565 

0.04 1 

233.9 

4.1178 

1678 

0 

1579 

1 1569 

0.047 

213.6 

3.7297 

1663 

0 

1578 

11574 

0.068 

218.9 

5.0287 

1662 

0 

1588 

20019 

0.090 

141.7 

4.9474 

1697 

0 

1636 

20056 

0.046 

225.7 

6.7501 

1687 

0 

1593 

20095 

0.057 

183.0 

4.3372 

1669 

0 

1604 

20096 

0.049 

182.7 

4.0470 

1682 

0 

1611 
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APFENDIX-3 


Sample Data set oi'planl (3) (Ibr details oTvaiiables see Appendix 1) 



O.OSd 
0.047 
0.043 
0.081 
0.057 
0.057 
0.008 
0.006 
0.142 
0.050 
0.088 
0.007 
0.055 
0.008 
0.103 
0.064 
0.064 
0.075 
0.076 
0.047 
0.057 
0.052 
0.046 
0.052 
0.043 
0.059 
0.049 
0.007 
0.04 I 
0.053 
0.051 
0.052 
0.0.53 
0.055 
0.038 
0.001 
0.082 
0.001 
0.043 
0.055 


155.67 

82.07 
128.00 

129.20 

118.53 

1.54.13 

1 13.40 

108.13 

137.07 

129.07 

147.53 

124.33 

129.07 

124.13 

140.67 
145.60 

121.33 
132.00 

157.27 

150.53 

125.13 

148.13 

129.07 

125.47 

134.53 

134.20 

180.40 

134.07 

147.47 

133.40 

127.27 
121.80 

147.40 
1 44. .3.3 

150.47 
I2(i.07 

1.54.40 

159.47 

140.27 
110.87 


3.3112 

3.7524 

5.5732 

4.0937 

6.6491 

4.2895 

3.7188 

3.9241 

3.7164 

3.4888 

4.3149 

5.6757 

3.721.5 

5.6294 

4.8219 

4.3692 

2.8124 

6.1096 

4.5839 

4.4429 

5.5283 

4.5990 

3.6186 

2.7817 

3.0713 

3.0031 

4.8102 

5.3248 

3.9.531 

3.4251 

5.1502 

5.2489 

4.7865 

5.0425 

4.4507 

0.3197 

4.0927 

3.8174 

2.9743 

4.5456 
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Sample Data set of plant (3) (continued .. ) 



Mmn 

HD 


0 22 

m 

Ore 2 

(kg) 

:i()(i;i 1 

3 i 6500 

0.516 

1 1.588 

0.477 

2.335 

0 

dOG.'f'l 

313900 

0.433 

1 0929 

0.122 

12.31 

0 

30();i7 

317100 

0.285 

10793 

0.277 

1920 

0 

30639 

327000 

0.384 

11015 

0.359 

1938 

0 

GOGdO 

327600 

0.204 

9522 

0.281 

1778 

0 

3()6dG 

326000 

0.317 

9499 

0.440 

2312 

0 

30647 

326400 

0.386 

10040 

0.377 

1701 

0 

30654 

326200 

0.468 

12837 

0.3.54 

1622 

0 

30(i65 

315500 

0.526 

13216 

0.636 

2056 

0 

30667 

324200 

0.464 

1 1246 

0.389 

19.36 

0 

30675 

313700 

0.364 

10558 

0.,546 

2213 

0 

30680 

328300 

0.237 

9463 

0.297 

1865 

0 

30682 

:V289()0 


9.521 

0.417 

1936 

0 

30683 

318800 

0.265 

10191 

0.375 

1862 

0 

30712 

316000 

0.295 

9632 

0..521 

2110 

0 

30713 

331000 

0.311 

9638 

0.377 

2184 

0 

30726 

321100 

0.857 

16584 

0.335 

1820 

0 

30753 

333200 

0.241 

10513 

0.276 

1980 

0 

30764 

324800 

0.406 

12953 

0.416 

2359 

0 

30767 

314100 

0.364 

10885 

0.359 

2258 

0 

30777 

331600 

0.294 

11549 

0.275 

1877 

0 

30806 

324300 

0.400 

12784 

0.362 

2222 

8 

30811 

330400 

0.368 

9428 

0.322 

1936 

0 

30812 

330400 

0.717 

14121 

0.341 

1882 

3 

30818 

334500 

0.602 

13253 

0.287 

2018 

0 

30821 

329300 

0.571 

12100 

0..349 

2013 

0 

30823 

331600 

0.495 

16919 

0.537 

2706 

0 

30844 

333700 

0.393 

14964 

0.298 

2011 

0 

30875 

332400 

0.476 

13403 1 

0.321 

2212 

0 

30876 

3 1 5700 

0.503 

1 16.55 

0.4,36 

2001 

0 

30888 

333400 

0,153 

16668 

0.247 

1909 

0 

30890 

314400 

0.353 

12483 

0.351 

1827 

0 

30897 

327900 

0.402 

13520 

0.400 

2211 

0 

30898 

333700 

0.351 

14162 

0.397 

2165 

0 

30910 

330500 

0.435 

13730 

0.441 

2347 

0 

30915 

329000 

0.213 

9490 

0.229 

1900 

0 

30921 

313500 

0.433 

11905 

0.642 

2316 

0 

30958 

328400 

0.411 

11041 

0.429 

2392 

0 

30992 

324700 

0.560 

11.589 

0.375 

2104 

0 

31013 

314500 

0.306 

9374 

0.374 

1663 

0 
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APPENDIX-4 


Sample Data set of plant (4) ' 


(for details of variables see Appendix 1) 


Heat No. 

0 22 

Tf 

Scrap 

Ore2 

HM 

lime 


(m')' 

(“C) 

(%) 

(kg) 

; (kg) 

(kg)' 

515991 

1380 

1661 

190000 

■602 

217510 

6310 

515992 

1150 

1650 

190300 

1597 

217420 

4650 

515993 

1150 

1648 

183300 

0 

224705 

8080 

515994 

2340 

1614 

181800 

2804 

223420 

15460 

515995 

1660 

1626 

187200 

979 

214305 

12370 

515996 

1330 

1639 

183100 

1329 

218635 

13550 

515997 

1310 

1659 

175300 

0 

214450 

21930 

515998 

1130 

1649 

169200 

0 

214600 

19700 

515999 

1420 

1658 

181200 

1006 

215600 

13260 

516000 

1550 

1629 

173100 

0 

212555 

19640 

516001 

1090 

1639 

175200 

0 

214900 

10380 

516002 

780 

1657 

189700 

1183 

222300 

8570 

516003 

1510 

1 637 

189400 

0 

222700 

10320 

516004 

1510 

1644 

192900 

1101 

226820 

13050 

516005 

1260 

1637 

192900 

1499 

228615 

4540 

516006 

1330 

1625 

183900 

1136 

216710 

7820 

516008 

120 

1654 

188100 

0 

220725 

9390 

516009 

1210 

1632 

188100 

1515 

218330 

11140 

516010 

1330 

1633 

187100 

1821 

221015 

7620 

516011 

1690 

1639 

189600 

2950 

222215 

9640 

516012 

1520 

1604 

189300 

2252 

222905 

8660 

516013 

1640 

1636 

187200 

2467 

219505 

7990 

516014 

1250 

1614 

188700 

1039 

221720 

9230 

' 516015 

0 

1634 

188400 

1809 

221920 

12660 

516016 

2120 

1617 

186600 

0 

223340 

16720 

516018 

1280 

1663 

183100 

3738 

221945 

4470 

516019 

1090 

1641 

185400 

3236 

219320 

5540 

516020 

1180 

1628 

182700 

2383 

220420 

7490 

516021 

1610 

1629 

190300 

3025 

223990 

6070 

516022 

1080 

1646 

181500 

2082 

219010 

6760 

516023 

1270 

1639 

179000 

2268 

217705 

7690 

516024 

910 

1643 

179900 

670 

217555 

4100 

516025 

1360 

1655 

181300 

1482 

218960 

4780 

516026 

1180 

1622 

183400 

1192 

222625 

5320 

516027 

1450 

1625 

183300 

972 

216735 

8090 

516028 

1620 

1629 

183500 

1677 

219780 

13070 

516029 

1680 

1622 

178600 

1683 

217260 

8490 

516030 

1540 

1641 

186000 

2049 

216210 

8190 

516031 

1680 

1607 

180000 

600 

222410 

17840 

516032 

1210 

1634 

188400 

0 

223610 

10340 


(continued ...) 


Sample Data set of' plant (4) (continued .. ) 




Si 

in HM 
% 

■ 

■ 


515991 

1584 

0.463 

0.044 

0.487 

190.3361 

515992 

1612 

0.169 

0.046 

0.377 

158.6791 

515993 

1584 

0.581 

0.031 

0.151 

153.5346 

515994 

1557 

0.585 

0.034 

0.764 

314.2064 

515995 

1550 

0.511 

0.044 

0.401 

232.3791 

515996 

1579 

0.553 

0.032 

0.250 

182.4959 

515997 

1565 

0.837 

0.031 

0.296 

183.2595 

515998 

1577 

0.861 

0.035 

0.217 

157.9683 

515999 

1583 

0.598 

0.032 

0.450 

197.5881 

516000 

1552 

0.779 

0.042 

0.474 

218.7669 

516001 

1580 

0.(557 

0.039 

0.273 

152.1638 

516002 

1612 

0.555 

0.049 

0.250 

105.2632 

516003 

1547 

0.586 

0.082 

0.709 

199.7178 

516004 

1567 

0.582 

0.051 

0.628 

165.3435 

516005 

1586 

0..336 

0.042 

0.462 

184.1170 

516006 

1571 

0.362 

0.0.56 

0.576 

180.5307 

516008 

1637 

0.499 

0.077 

0.077 

228.1574 

516009 

1584 

0.327 

0.070 

0.499 

224.1407 

516010 

1583 

0.280 

0.040 

0.537 

169.1322 

516011 

1583 

0.473 

0.060 

0.589 

173.0158 

516012 

1574 

0.594 

0.067 

0.462 

149.0972 

516013 

1579 

0.482 

0.033 

0.474 

215.6346 

516014 

1573 

0.528 

0.038 

0.250 

147.9384 

516015 

1561 

0.325 

0.037 

0.537 

175.0075 

516016 

1523 

0.654 

0.092 

0.984 

125.4855 

OC 

o 

1628 

0.388 

0.041 

0.589 

186.3354 

516019 

1622 

0.399 

0.054 

0.462 

159.0118 

516020 

1602 

0.613 

0.073 

0.589 

200.7059 

516021 

1580 

0.478 

0.046 

0.550 

221.1302 

516022 

1608 

0.462 

0.041 

0.365 

231.9801 

516023 

1601 

0.592 

0.0.54 

0.499 

213.6811 

516024 

1601 

0.244 

0.045 

0.296 

162.3362 

516025 

1598 

0.281 

0.034 

0.450 

234.8309 

516026 

1572 

0.267 

0.034 

0.273 

205.6386 

516027 

1661 

0.493 

0.038 

0.537 

237.0572 

516028 

1573 

0.552 

0.059 

0.628 

231.1207 

516029 

1559 

0.600 

0.050 

0.709 

207.5133 

516030 

1580 

0.601 

0.043 

0.589 

205.7707 

516031 

1551 

0.626 

0.()65 

0.576 

156.6082 

5 I ()032 

1548 

0.399 

0.036 

0.330 

226.8765 


(continued ...) 
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12^1 


Sample Data set of plant (continued ...) 


) 





■ 

cloloinilc. 

AOj 
my ton 

'(Cr-Cst) 

(l/Ci-l/Csi.) 

515991 

2820 

7.8189 

3.3474 

9130 

6.3445 

477 

-0.0021 

515992 

2370 

7.4336 

6.7474 

7020 

5.2893 

367 

-0.0027 

515993 

4330 

4.7984 

3.5406 

12410 

5.1178 

141 

-0.0066 

515994 

2120 

4.6735 

6.7837 

17580 

10.4735 

754 

-0.0013 

515995 

3580 

7.8326 

6.0346 

15950 

7.7460 

391 

-0.0025 

515996 

3760 

5.7943 

6.2450 

17310 


240 

-0.0040 

515997 

4260 

4.8966 

6.9749 

26190 


286 

-0.0034 

515998 

4070 

4.0190 

6.3106 

23770 


207 

-0.0046 

515999 

3730 

5.8452 

5.7107 

16990 


440 

-0.0022 

516000 

4130 

4.7950 

6.7969 

23770 


464 

-0.0021 

516001 

4100 

4.8000 

4.2083 

14480 


263 

-0.0037 

516002 

4200 

6.0997 

3.7986 

12770 


240 

-0.0040 

516003 

4060 

6.0470 

5.4245 

17110 


699 

-0.0014 

516004 

3140 

5.4801 

3.2688 

7680 


618 

-0.0016 


3150 


5.4818 

10970 


452 

-0.0022 


4440 


6.7878 

12060 


566 

-0.0017 


4000 


5.0163 

13640 


489 

-0.0020 

516009 

4170 


4.1324 

12160 


527 

-0.0019 

516010 

4290 


4.3232 

13520 


579 

-0.0017 

516011 

3120 

5.1577 

2.9363 

7590 


452 

-0.0022 

516012 

3410 

5.4052 

3.4949 

8950 

7.4714 

464 

-0.0021 

516013 

2940 

5.8916 

3.1141 

9010 

5.6377 

240 

-0.0040 

516014 

3550 

5.1709 

3.7621 

10310 

9.4923 

974 

-0.0010 

516015 

3270 

4.9311 

3.3866 

10960 

5.7672 

579 

-0.0017 

516016 

2390 

4.9972 

4.3588 

6490 

4.9699 

452 

-0.0022 

516018 

2350 

5.0361 

4.3786 

7130 

5.3534 

579 

-0.0017 

516019 

2380 

4.9037 

5.0700 

7700 

7.1878 

540 

-0.0018 

516020 

3830 

5.3754 

4.1778 

11920 

4.9313 

355 

-0.0027 

516021 

1530 

5.2730 

6.0215 

14600 

5.8336 

489 

-0.0020 

516022 

3890 

4.9066 

3.6973 

12380 

4.1829 

286 

-0.0034 

516023 

3860 

6.6667 

3.4190 

12050 

6.2112 

440 

-0.0022 

516024 

3520 

5.6071 

6.4191 

13860 

5.3004 

263 

-0.0037 

516025 

3390 

5.7066 

7.2260 

18320 

6.6902 

527 

-0.0019 

516026 

2980 

6.2561 

7.5838 

16760 

7.3710 

618 

-0.0016 

516027 

3620 

6.2145 

9.4463 

24050 

7.7327 

699 

-0.0014 

516028 

3900 

4.9429 

2.6490 

9610 

7.1227 

579 

-0.0017 

516029 

2950 

6.7071 

2.7488 

9200 

7.5536 

566 

-0.0017 

516030 

2880 

6.7082 

2.8634 

8640 

5.4112 

320 

-0.0030 

516031 

3020 

6.6585 

3.1760 

10330 

7.8277 

440 

-0.0022 

516032 

3800 

5.6585 

3.1862 

10250 

6.8546 

355 

-0.0027 












